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ABSTRACT 

We present a Nobeyama 45 m Radio Telescope map and Australia Telescope Com- 
pact Array pointed observations of N 2 H + 1-0 emission towards the clustered, low mass 
star forming Oph B Core within the Ophiuchus molecular cloud. We compare these 
data with previously published results of high resolution NH3 (1,1) and (2,2) observa- 
tions in Oph B. We use 3D clumpfind to identify emission features in the single-dish 
N2H + map, and find that the N 2 H + 'clumps' match well similar features previously 
identified in NH3 (1,1) emission, but are frequently offset to clumps identified at similar 
resolution in 850 /im continuum emission. Wide line widths in the Oph B2 sub-Core 
indicate non-thermal motions dominate the Core kinematics, and remain transonic at 
densities n ~ 3 x 10 5 cm -3 with large scatter and no trend with JVfHk). In contrast, 
non-thermal motions in Oph Bl and B3 are subsonic with little variation, but also show 
no trend with H2 column density. Over all Oph B, non-thermal N 2 H + line widths are 
substantially narrower than those traced by NH3, making it unlikely NH3 and N 2 H + 
trace the same material, but the v LSR of both species agree well. We find evidence for 
accretion in Oph Bl from the surrounding ambient gas. The NH3 / N 2 H + abundance 
ratio is larger towards starless Oph Bl than towards protostellar Oph B2, similar to 
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recent observational results in other star-forming regions. The interferometer observa- 
tions reveal small-scale structure in N2H" 1 " 1-0 emission, which are again offset from 
continuum emission. No interferometric N2H + emission peaks were found to be coinci- 
dent with continuum clumps. In particular, the ~ 1 M B2-MM8 clump is associated 
with a N2H + emission minimum and surrounded by a broken ring-like N2H + emission 
structure, suggestive of N2H 4 " depletion. We find a strong general trend of decreasing 
N2H + abundance with increasing iV(H2) in Oph B which matches that found for NH3. 

Subject headings: ISM: molecules - stars: formation - ISM: kinematics and dynamics - 
ISM: structure - radio lines: ISM 



1. Introduction 



Stars form out of the gravitational collapse of centrally condensed cores of dense molecular 
gas. Recent years have seen leaps forward in our understanding of the structure and evolution 
of isolated, star for ming cores. Most star formation, however, occurs in clustered environments 
(jLada Ladall2003l ). These regions are more complex, with complicated observed geometries, and 
contain c ores which tend to have high er densities and more compact sizes than those found in 
isolation (["Ward-Thompson et al.ll2007al ). It is likely that due to these differences the evolution of 
filaments and cores in clustered regions proceeds differently than in the isolated cases. Character- 
izing the physical and chemical structures of these more complicated regions are thus the first steps 
towards a better understanding of the process of clustered star formation. 

It is now clear that molecular cores become extremely chemically differentiated, as many 
molecules commonly used for tracing molecular gas, such as CO, bec ome severely depleted in the 
innermost core regions throug h adsorption onto d ust grains (see, e.g.. lDi Francesco et al.ll2007l . for 
a review). In a recent paper rtFriesen et al.ll2009l . hereafter Paper I), we studied the dense gas in 
several cluster forming Coreqj in the Ophiuchus molecular cloud (Oph B, C and F) through high 
resolution observations of NH3 (1,1) and (2,2) emission made at the Green Bank Telescope (GBT), 
the Australia Telescope Compact A rray ( ATC A) and the Very Large Array (VLA) . The Ophiuchu s 
molecular cloud (~ 120 pc distant; lLoinard et al.ll2008l : lLombardi et al.ll2008l : LKnude Hog)ll998l ) 
is our closest example of ongoing, clustered star formation. We found that the Ophiuchus Cores 
presented physical characteristics similar to those found in studies of both isolated and clustered 
environments. In Oph C, for example, gas motions are subsonic (mean a v = 0.16 kms -1 ), and 



1 In this paper, we describe Oph A, B, C, etc ., as 'C ores' since this is how these features were named in DCO + 
observations of the L1688 region bv lLoren etaD (|l990h . Since then, higher-resolution data such as those described 
in this paper have revealed substructure in these features that could be themselves precursors to stars, i.e., cores. 
To avoid confusion, we refer the larger features in Oph as 'Cores' and Core substructure identified by CLUMPFIND as 
'clumps'. 
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decrease in magnitude along with the kinetic gas temperature (Tr-) towards the ther mal dust con- 



tinuu m emission peak in a manner reminiscent of findings in the isolated cores L1544 (jCrapsi et al 



20071 ); etc. In contrast to Oph C, NH3 line widths in Oph B are dominated by transonic non- 
thermal motions and the gas temperatures are warmer than typically found in isolated regions 
((Tr) = 15 K) and nearly constant across the Core. No contrast in any parameters [such as Tr-, 
ratio of non-thermal to thermal line width a NT / c s , or NH3 column density iV(NH3)] was found 
between the two main components of Oph B, labelled Bl and B2, despite the presence of several 
embedded Class I protostars in B2 while Bl appears starless. Additionally, on small scales (~ 15", 
or 1800 AU at 120 pc), significant offsets were found between the peaks of integrated NH3 (1,1) 
intensity and those of submillimete r continuum emission , as well as between individual 'clumps' 
identified through 3D clumpfind ([Williams et al.l Il994l ) in the NH3 data cube and those found 
through 2D clumpfind in 850 p continuum emission. Finally, evidence was found for a decreas- 
ing fractional NH3 abundance with increasing H2 column density (as traced by 850 /um continuum 
emission), suggesting that NH3 (1,1) may not be tracing the densest gas in the Oph Cores. 

The diazenyl ium ion, N 2 H + , has been shown to be a preferential tracer of quiescent dense gas in 



2004 



Di Francesco et al 



mole cular clouds (IWomack et alll992l : lCaselli et al.ll2002al : lTafalla et al.ll2002l . 
2004). The 1-0 rotational transition has a critical density n CT ~ 2 x 10 5 cm" 3 , ~ 10 — lOOx greater 
than that of NH3 (1,1), and is therefore possibly better suited to probing gas properties at the high 
densities expected in clustered star forming Cores. Based on both observations and chemical mod- 
els of starless cores, N2H + also appears resilient to depletion at the high densities (n > 10 5 cm" 3 ) 



Tafalla et al. 


2004; 


Bergin et al. 


2002) 



e.g., 

physical conditions of dense cores and clumps. 



Recently, lAndre et al.1 ([20071 ) observed the Oph A, B, C, E and F Cores in N 2 H + 1-0 emis- 
sion at 26" angular resolution using the IRAM 30 m telescope. They foun d that the gas trac ed 
by N2H + emission near clumps identified in millimeter continuum emission (jMotte et al.lll998l ) is 
characterized, on average, with relatively narrow line widths showing small non-thermal velocity 
dispersions, in contrast to the large non-thermal motions traced by NH3 emission in Paper I. The 
authors derived virial masses from the N2H + emission, and found they generally agreed within a 
factor ~ 2 with the masses derived from dust emission, suggesting that the clumps are gravitation- 
ally bound and prestellar. The relative motions of the clumps are also small and subvirial, with a 
crossing time larger than the expected lifetime of the objects, such that clump-clump interactions 
are not expected to impact their evolution. These results suggest that a dynamic picture of clump 
evolution involving competitive accretion at the prestellar stage does not accurately describe the 
star formation process in central Ophiuchus. 

In this work, we discuss the results of higher resolution (18", or ~ 2200 AU at 120 pc) observa- 
tions of N2H + 1-0 in the Ophiuchus B Core that reveal the distribution, kinematics and abundance 
pattern of the Core and associated embedded clumps on smaller physical scales. We additionally 
examine N2H" 1 " 1-0 structure at even higher resolution (8" x 5"), through observations made with 
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the ATCA, towards five locations in Oph B2 where small-scale structure was found in NH3 (Paper 
I) and in previously unpublished Berkeley-Illinois-Maryland Association (BIMA) N2H + 1-0 obser- 
vations. We compare these data with our recently published analysis of NH3 emission in Oph B, in 
particular focusing on the kinematics and relative abundances of the two species to compare with 
physical and chemical models of dense core formation and collapse. 

We discuss the observations and calibration and §2. In §3, we present the data, and discuss 
the results of the hyperfine line structure fitting procedure and derivations of the column density, 
iV(N2H + ), and fractional N2H + abundance, X(N2H + ), in §4. In §5, we discuss general trends in 
the data, and compare these results with those found in Paper I and with studies of dense cores in 
isolated environments. We summarize our findings in §6. 



2. Observations and Data Reduction 
2.1. Nobeyama 45 m Radio Telescope 

Single-dish mapping observations of Oph B in N2H + 1-0 (rest frequency 93.1 74 GHz) were 



performed May 10 - 15, 2007 using the 25-BEam Array Receiver System (BEARS, ISunada et al 



2000; lYamaguchi et al.ll2000l ) at the Nobeyama 45 m Ra dio Telescope. BE ARS is a 5 x 5 SIS het 



erodyne receiver array. The autocorrelator spectrometer (jSorai et al.ll2000l ) with 8 MHz bandwidth 
and 1024 channels was used as the backend, giving a spectral resolution of 7.8 kHz or 0.025 km s -1 
at 93 GHz. The final spectral resolution was binned to 0.05 kms -1 . The half-power beamwidth was 
~ 18" and the main-beam efficiency at 93 GHz is estimated to be tjmb = 0.505, interpolated from 
observatory measurements at 86 GHz and 100 GHz, with an uncertainty of ~ 8 — 10 %. The beams 
are separated by 41.1" on the sky. 

Observations were performed in position-switching OTF mode, alternating scanning in Right 
Ascension (R.A.) and Declination (deck) to ensure minimal striping artifacts. The array was 
rotated by an angle of 7° such that the spacing of individual detectors was 5" in R.A. or deck, 
and scans of the array were spaced by 25". A ten minute integration on the off-source position 
(+5° in R.A.) was performed to ensure no emission was present. The weather was only moderately 
variable over the observations, and most data were taken with T sys ~ 300 — 450 K (DSB). Telescope 
pointing was checked every 90 minutes on nearby SiO maser sources AF Sco and RS Lib. Typical 
pointing offsets were 1 — 3". 

Calibration of BEARS was done by observing a bright N2H + target first with the S100 single 
sideband receiver, and second with each individual BEARS beam. These data provided a flux 
scaling factor for each beam to correct for the telescope gain and DSB ratio. The starless core 
L183 was observed with S100 and the central BEARS beam, but the line emission was not strong 
enough to use as a calibrator for each BEARS beam. The Oph A core SMI was then observed with 
each BEARS beam and the scaling factors for each beam were determined using both sources. The 
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average scaling factor was 1.35, with a range from 1.07 to 1.66. 

The scaling factors were applied to the data and a linear baseline was fit to and subtracted 
from the spectra. The R.A. and decl. scans were grid ded separately and then c ombined using a 
basket- weave method to minimize baseline drift effects (jEmerson &: Graevd ll988). 



2.2. BIMA 

Observations of N 2 H+ 1-0 across the Oph B2 Core were made with the BIMA array located at 
Hat Creek, CA, U.S.A. Four tracks were observed on 2001 March 18, 21, and 26 and 2001 May 23, 
when the array was in its C-configuration. Observations of Oph B2 consisted of a serial sequence of 
12 overlapping positions with centers spaced by 1.06', to obtain a Nyquist-sampled mosaic across 
Oph B2 at 93.2 GHz. The correlator was configured to provide 12.5 MHz of bandwidth in one 
correlator window at 93.2 GHz, providing 512 channels of 0.04 km s™ 1 width. The radio source 
3C279 was observed at the start of each track to provide passband and flux calibration. Neptune 
was observed at the end of each track to provide flux calibration. The radio source 1733-130 was 
observed between every set of Oph B2 observations to provide phase calibration. Values of T sys 
during the tracks were typically 250-500 K for all antennas, except antenna 2 where T sys values 
were a factor of ~ 1.5 higher. 



The Oph B2 BIMA data were reduced using standard tasks in the MIRIAD package (jSault et al 



19951 ). To improve sensitivity to diffuse emission and produce a rounder beam, the visibilities were 
tapered with a Gaussian 10" FWHM in R.A. and 0.1" FWHM in decl. For the N 2 H+ data, the 
final beam size achieved was 19.3" x 10.6" FHWM (P.A. -12.4°) and the resulting sensitivity (1 
a rms) per channel was 0.72 Jybeam -1 , or equivalently 0.49 K. Only faint, diffuse N 2 H + emission 
was detected towards locations in eastern Oph B2, with even weaker emission spots towards central 
Oph B2. The faint N 2 H + emission was not detected significantly enough to fit the hyperfine line 
components, but the locations of significant integrated intensity provided targets for later, more 
sensitive data of Oph B2 with the ATCA (see below). 



2.3. ATCA 

A total of 11 pointings were observed in N 2 H + 1-0 emission at the ATCA towards five targets 
in July 2006 in the interferometer's most compact (H75) configuration. The BIMA data were used 
as a guide for the pointings. Figure shows the positions of all 11 pointings over an 850 /im 
continuum map of Oph B. Table Q] lists the positions of the targets, the rms noise achieved in 
Jybeam -1 , and any associated NH3 clumps identified in Paper I, continuum objects, or embedded 
Young Stellar Objects (YSOs). Two of the targets were imaged using mosaics of three Nyquist- 
spaced pointings each. Three pointings towards the Oph B2 continuum emission peak overlapped, 
but were not Nyquist-spaced, while the remaining two targets were observed with a single pointing 
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each. An 8 MHz band width with 512 channels was used for the observations, providing a spectral 
resolution of 15.6 kHz or 0.05 kms -1 at 93.2 GHz while allowing all hyperfine components of the 
transition to fit within the band. Baselines ranged from 31m to 89 m for five antennas (the sixth 
antenna is located 6 km from the main group and was not used for these observations due to poor 
phase stability), corresponding to 10.3 kA - 29.7kA. The primary beam (field-of-view) was 34" in 
diameter. 

Observations were made in frequency-switching mode. Observations cycled through each in- 
dividual pointing between phase calibrator observations to maximize u, v coverage and minimize 
phase errors. The phase calibrator, 1622-297, was observed every 20 minutes, and was also used to 
check pointing every hour. Flux and bandpass calibration observations were performed each shift 
on Mars or Uranus. 



The ATCA data were reduced using the Miriad data reduction package (jSault et al.lll995l ). The 
data were first flagged to remove target observations unbracketed by phase calibrator measurements, 
data with poor phase stability or anomalous amplitude measurements. The bandpass, gains and 
phase calibrations were applied, and the data were then deconvolved (or jointly deconvolved, in the 
case of the mosaiced pointings). First, the data were transformed from the u, v plane into the image 
plane using the task INVERT and natural weighting to maximize sensitivity. The data were then 
deconvolved and restored using the tasks MOSSDI and RESTOR to re move the beam re sponse 



from the image. MOSSDI uses a Steer-Dewdney-Ito CLEAN algorithm (jSteer et al.lll984l ). The 
cleaning limit was set at twice the rms noise level in the beam overlap region for each object. Clean 
boxes were used to avoid cleaning noise in the outer regions with less beam overlap. The final 
synthesized beam had a FWHM of 8" x 5". 

We next convolved the ATCA data only to a final beam FWHM of 18" to match the Nobeyama 
data and estimate the amount of flux recovered by the interferometer. With the exception of 
the single pointing towards the NH3 clump B2-A7, we found the maximum peak intensities in 
the convolved maps are only ~ 15 — 20 % of that in the single-dish map, suggesting that the 
interferometer data alone are only sensitive to ~ 15 — 20 % of the total N2H" 1 " emission. At the 
B2-A7 location, the single detected interferometer emission peak is > 50 % of that seen by the 
single-dish telescope. 



Results 



3.1. Single-Dish Data 



In Figure [T^,, we show 850 /im continuum emission in Oph B at 15" FWHM resolution, first 
mapped with the Submillimet re Common Use r Bolom eter Array (SCUBA) at the James Clerk 
Maxwell Telescope (JCMT) by I Johnstone et al.1 (l2000br) and recen t ly re- reduced and combined with 
all other SCUBA archive data in the region by iJorgensen et al.1 (|2008l ). following the description 
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iii 



Kirk et al.l ([2006) . The Oph B Core can be separated into three distinct sub-Cores, Bl, B2, 



and B3, which are labelled on the Figure. We also show the locations of four C lass I protostars 



(Vrba et al. 


1975: 


Elias 


1978) 



Figure [Tb shows the N2H + 1-0 intensity in the Oph B Core at 18" FWHM resolution, inte- 
grated over only those velocity channels containing emission from the hyperfine components. On 
18" scales, the integrated N2H + intensity within Oph Bl and B2 is remarkably smooth and shows 
several maxima on scales significantly larger than the beam. On average, the integrated N2H + 
intensity closely follows the continuum emission in the Core, with a few exceptions. First, we find 
significant N2H + emission in southeast Oph B2, where little continuum emission is seen. In fact, 
the entire eastern half of B2 contains strong N2H + emission, while the continuum emission becomes 
significantly fainter to the east of the central B2 850 /im emission peak. Second, the integrated 
N2H + intensity maximum in B2 is offset from the continuum maximum by ~ 23", although strong 
N2H + emission is also found at the continuum peak. Third, the continuum and N2H + integrated 
intensity peaks in Oph Bl are also offset. Fourth, the filament connecting Oph Bl and B2 is more 
prominent in N2H + than in continuum emission. The Oph B3 Core is undetected in 850 fim con- 
tinuum. Although it is a clear feature in the N2H + data, it is faint in N2H + integrated intensity 
due to the extremely narrow line widths associated with the Core (discussed further in 



3.1.1. Identifying discrete N^IP" 'clumps' 
We next identify discrete emission 'clumps' in the Nobeyama N2H + 1 -0 emission data cub e 



using the 3D version of the automated structure-finding routine clumpfind (jWilliams et al.lll994l ). 
CLUMPFIND uses specified brightness contour intervals to search through the data cube for distinct 
emission features identified by closed contours in position and velocity. The size of the structures 
found are determined by including adjacent pixels down to an intensity threshold, or until the 
outer edges of two separate structures meet. The clump location is defined as the emission peak 
within the core. The identification of emission clumps in the data allows us to discuss the absolute 
locations of N2H + emission within Oph B with respect to the locations of objects identified in 
continuum emission, as well as emission features in other molecular line tracers. We caution that 
the identification of N2H + 'clumps' in Oph B does not imply any association with underlying 
physical structure in the core. 

The rms noise level in the data cube is essentially constant (rms ~ 0.08 K in Tt units) over 
most of the region observed, but emission in Oph Bl extends near the edge of the area of best OTF 
mapping overlap and consequently has slightly lower sensitivity. We therefore created a signal-to- 
noise data cube by dividing the spectrum at each pixel by the rms noise in the off-line spectral 
channels. We then ran clumpfind specifying brightness contours of 4x the S/N ratio with an 
intensity threshold of 3x the data rms. The strongest hyperfine components of the N2H + 1-0 line 
overlap significantly over much of the region observed, so we ran the algorithm on the isolated 
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F\F — > F[F' = 01 —* 12 component only. This strategy allowed the clear separation of multiple 
velocity components along the line of sight, but reduced the signal-to- noise ratio of the detections 
since the isolated component can be substantially fainter (up to a factor of ~ 2.3) than the central 
hyperfine emission. We therefore performed several checks to ensure identified N2H + objects were 
legitimate. We removed any identified clumps with areas less than the Nobeyama beam, clumps 
with a peak location at the edge of the mapped region, and clumps with reported line widths less 
than twice the velocity resolution of the observations (0.05 kms -1 ). These limits culled most of 
the spurious objects, and we finally removed several additional objects through visual inspection. 

We found a total of 17 separate N2H + clumps in Oph B from the Nobeyama data. Table 
[2] lists the locations and FWHM in AU of the clumps, as well as the peak temperatures (T|) of 
the isolated and main hyperfine N2H + components, and the signal-to-noise ratio of the isolated 
hyperfine component. We plot the positions of the clumps, i.e., the positions of peak brightness 
temperature of the isolated component, in Figure QJ). The centroid velocities and line widths of the 
N2H + clumps are also returned by clumpfind, but we determined these values instead through 
spectral line fitting (described in further detail in §3). The line fitting algorithm provided more 
accurate results given the hyperfine structure of the N2H + 1-0 line. 

We show in Figure [2] the N2H + clumps as ellipses overlaid on the emission at the velocity 
channels that most closely match the clumps' centroid velocities determined through spectral line 
fitting in §3. The ellipse axes are the returned clump FWHMs in R.A. and deck The 850 /.mi 
continuum emission contours are also shown in Figure [2j In many cases, an offset between the 
N2H + clumps and continuum emission is clear. 

We first compared th e locations of the N2H" 1 " clumps with the peak locations of 850 |im contin- 



uum clumps identified bv lJorgensen et al.l (|2008l ) in the SCUBA data. Figured shows the angular 
separation between N2H + clumps and the peak location of the nearest 850 /xm continuum clump. 
Few continuum clumps lie within a radius of 18" in angular distance from an N2H + clump, sug- 
gesting that many of the N2H + clumps thus identified are not correlated with continuum objects. 
This result is similar to our findings in Paper I, where a similar analysis of NH3 emission in Oph B 
revealed significant offsets between locations of NH3 clumps and continuum clumps. Interestingly, 
the mean minimum separation between N2H + and continuum clumps in Oph Bl (~ 20") is half 
that in Oph B2 (~ 40"), suggesting there is a better correlation between N2H + and continuum 
emission in Oph Bl. Furthermore, Figure [3b shows the angular separation between N2H + clumps 
and nearest NH3 clumps from Paper I, and shows a clear correlation between the two sets of ob- 
jects. In the majority of cases, N2H + and NH3 clumps are coincident within the resolution of the 
N 2 H+ data. 

Other searches for submillimeter continuum emission clumps in Oph B using different clump- 
finding algorithms produce different num bers of objects within each sub-Core. For example, the 



clump catalogue by ISimpson et akl (|2008l ) was created fro m a combination of all 850 u m SC UBA 



data in Ophiuchus, following the technique described by iNutter &: Ward-Thompson! (|2007l ). and 
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cont ains a larger n umber of objects within the Oph B Core. Comparison of the N2 H" 1 " clumps with 



the lSimpson et al J continuum clumps shows a better correspondence than with the I Jgrgensen et al 



catalogue, but the median offset between N2H + and continuum clump locations is still greater than 
the 18" FWHM beam of the N2H + observations. In the following, we restrict our analysis to the 
J0rgensen et al.l catalogue as we use the same clump identification technique, but note that offsets 
are also seen between the N2H + clumps presented here and other continuum clump catalogues. 

We furt her compared the lo cations of the N2H + clumps thus determined with N2H" 1 " clumps 
identified by lAndre et al.l (|20071 ) in N2H + 1-0 emission obser yed with the IRAM 30 m telescope 
(~ 26" half power beam width). We find good agreement with I Andre et al.l for most of the objects 
despite their use of a spatial filtering scheme rather than the CLUMPFIND algorithm to identify 
clumps. 



3.2. ATCA Data 



In Figures Hh, EK, Eti and[3 we show the integrated N2H + 1-0 intensity towards four of the five 
regions observed in Oph B2 with the ATCA. In the fifth, most western region observed (see Figure 
[Ik.), we found no significant interferometer emission above the rms noise. In all plots, the data were 
integrated over only those velocity channels containing emission from the hyperfine components. 
N2H + , NH3 and continuum clumps mentioned in the text are labelled on the Figures. 

Figure Hk shows the combined results of three ATCA pointings towards the 850 /i m continuum 



peak in Oph B2, with the locations of th e asso ciated continuum clump B2-MM8 (jMotte et al 



19981 . also 16273-24271 in I.Torgensen et al 



200. 



the nearby N2H + clump B2-N5 identified in 



3.1, and the Class I protostars Elias 32 ( Eliasl 1978 . also VSSG17; Vrba et al. 1975 ) and Elias 33 



(also VSSG18) identified (but note Elias 33 is outside the ATCA primary beam). The single-dish 
integrated intensity peak, coincident with B2-N5, is offset to the continuum data and B2-MM8 by 
~ 23" (see Figure [TJa). Two significant integrated intensity peaks are seen in the ATCA data. The 
first is coincident with B2-N5 (and consequently offset from B2-MM8), while the second is located 
between the two protostars where little continuum emission is seen. Some emission is also found 
to the east and west of B2-MM8. Negative emission features are coincident with the continuum 
clump peak, and additionally between the two strong N2H + integrated intensity features. The 
integrated intensity of the data forms a broken ring-li ke structure around B2-MM8. The distribution 
is reminiscent of N2H + 1-0 emission in B68, where iBergin et al.1 (|2002l ) found strong evidence of 
N2H + depletion through radiative tr ansfer modelling. Eviden ce for N2H + depletion on small scales 
has also been found in IRAM 04191 (jBelloche & Andrei I2OO 4 ). 



In Figure [5^, we plot the integrated intensity towards th e single point i ng in eastern B2, which 
includes the N 2 H+ clump B2-N6 and the NH 3 clump B2-A7. llVIotte et aP (|l998l ) i dentify a contin 



uum cl um p here (B2-MM15 , unre solved at 11"), but no object was identified by I Johnstone et al 



(|2000bl ) or |j0rgensen et al.1 (J2008J) in 850 fim continuum emission. In the interferometer data, a 
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single integrated intensity peak is seen slightly offset (~ 6", significantly less than the Nobeyama 
beam FWHM of 18") NNW from B2-N6. The emission is extended along the beam's long axis, 
but fitting a 2D Gaussian to the data indicates the emission is also unresolved by the 8"x 5" 
FWHM beam (~ 960 AU x 600 AU at 120 pc). The peak integrated intensity of the interferometer 
data is greater relative to the single-dish data in this location by a factor ~ 2 compared with the 
other interferometer observation locations, suggesting that a larger fraction of N2H + emission is 
associated with the small-scale peak rather than the extended N2H + emission. 

Figured shows the results of the south-eastern three pointing mosaic. Strong NH3 emission, 
including the NH3 clump B2-A9, was found in this region in Paper I. The only continuum emission 
in the mosaiced area is coincident with the Class I protostar Elias 33. Two overlapping N2H + 
clumps identified in §3.1.1 are shown, with clump peaks separated by ~ 8". We find two N2H + 
integrated intensity peaks in the region, separated by ~ 30", of which one is coincident with B2- 
N8. Two regions of negative integrated intensity are found within the primary beam contours. The 
most negative feature is coincident with Elias 33, while the second is near B2-N7. 

Figure [7] shows the combined results of the three-pointing mosaic in north-eastern Oph B2. 
One N2H + clump, B2-N9, is identified in the mosaiced area, which is offset from the continuum 
clump B2-MM16/ 16273-24262 by ~ 30". In the interferometer data, no significant emission is 
associated with either the N2H + or continuum clumps. One significant feature is visible in the 
Figure which is coincident with the western elongation (away from the continuum contours) of the 
single-dish N2H + integrated intensity. Multiple negative features are present. 

In general, the small scale structure visible in the integrated intensity maps does not correlate 
well with objects identified in continuum emission, and is also sometimes offset from clumps iden- 
tified on larger scales from the single-dish N2H + emission. The offsets between the single-dish and 
interferometric N2H + emission peaks tend to be small, however, and are likely due to the increase 
in spatial resolution afforded by the interferometer, which is better able to locate precisely the 
N2H + emission peak for a given object while resolving out large scale structure. No small-scale 
emission was found associated with embedded protostars. 



4. N2H + Line Analysis 



The N2H + 1-0 rotational transition has seven hyperfine components (ICaselli et al.lll995l ). We 
fit each observed N2H + 1-0 spectrum with a custom Gaussian hyperfine structure fitting routine 
written in IDL, described in detail in Paper I, assuming local thermodynamic equilibrium (LTE) and 
equal excitation temperatures T ex for each hyperfine component. The e xpected emission fre quencies 
and emitted line fractions for the hyperfine components were taken from lPagani et al.1 (120091 ). We fit 
the spectrum at each pixel and obtain maps of the line-of-sight local standard of rest (LSR) velocity 
(^lsr) °f the emission, the line FWHM (Av), the total line opacity (r), and the line brightness 
(Tmb)- We additionally calculate T ex assuming the line emission fills the beam (the determined 
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T ex will be a lower limit if the observed emission does not entirely fill the beam). 

The line widths determined by the hyperfine structure fitting routine are artificially broad- 
ened by the velocity resolution of the observations (0.05 kms -1 ). We therefore subtracted in 
quadrature the resolution width, Av res , from the observed line width, Av b s , such that Avu ne = 



The uncertainties reported in the returned parameters are those determined by the fitting 
routine, and do not take the calibration uncertainty (~ 8 — 10 % for the Nobeyama data) into 
account. The calibration uncertainty does not affect the uncertainties returned for v LSR or Av. The 
excitation temperature, N2H + column densities, fractional N2H + abundances and volume densities 
discussed further below, however, are dependent on the amplitude of the line emission, and are thus 
affected by the absolute calibration uncertainty. 

Over most of the map, the spectra could be fit with a single set of hyperfine components, but 
several small regions showed clear double-peaked line profiles in the isolated N2H + 1-0 hyperfine 
component, including the N2H + clumps Bl-Nl, B1-N3 and B1-N4. To understand better this 
emission, we fit the spectra with the above Gaussian hyperfine structure fitting routine using two 
separate velocity components. We show the spectra at the three N2H 4 " clump peak locations in 
Figure El overlaid with both single and double velocity component fits. In all cases, we are able to 
determine well the v LSR and Av of both components, but the fainter component is typically very 
optically thin and we are unable to constrain simultaneously both t and T ex . Although the angular 
separation of B2-N7 and B2-N8 clumps is small, visual inspection of the spectra at the clump peaks 
did not show a separate velocity component above the rms noise. We discuss the nature of the 
double-peaked objects further in §5.3.21 In the following sections, we discuss the results of the 
single velocity component fits only, unless otherwise stated. 



Table [3] lists the mean, rms, minimum and maximum values of v LSR , Av, r and T ex found 
in Oph B. The results across Oph B are shown in Figured! where we plot only data points with 
S/N ^ 2.5 in the isolated F 1 F -> F[F' = 01 ^ 12 component of the N 2 H+ 1-0 line. In the 
following sections, we describe in detail the results of the line fitting and examine the resulting 
line centres and widths, as well as T ex and t. Table 1^1 lists the determined parameters at the peak 
location of each N2H + clump. 



We show in Figures [9^ andEb, respectively, the fitted values of v LSR and Av of N2H + 1-0 across 
Oph B. Relatively little variation in line-of-sight velocity is seen across Oph Bl and B2, with no 




Avj, es . In the following, we simply use Av = Avu ne for clarity. 



4.1. Single-Dish Results 



4-. 1.1. Line Centroids and Widths 
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evidence of a discontinuity in velocity separating the two sub-Cores. A gradient of increasing v LSR 
from west to east is seen across Oph B, with the largest v LSR values found in the eastern half of B2, 
to a maximum v hSR = AAA kms -1 . In Oph B3, the mean v LSR = 3.06 kms -1 , significantly lower 
than in Bl ((v hSR ) = 3.97 kms -1 ) or B2 ((v LSR ) = 4.05 kms -1 ). The rms variation in v hSR across 
all Oph B is only 0.23 kms -1 , which is dominated by the spread in v LSR in B2 (rms = 0.20 kms -1 , 
compared with 0.12 kms -1 and 0.14 kms -1 in Bl and B3, respectively). 

We find greater variation between Oph Bl and B2 in N2H + 1-0 line width than in v LSR . Oph 
Bl is generally characterized by narrow lines, with a mean Av = 0.42 kms -1 and an rms variation 
&Av of only 0.08 kms -1 . Oph B2 contains regions of both small and large line widths, with a larger 
mean Av = 0.68 kms -1 and a^v = 0.21 kms -1 than found in Bl. Both Bl and B2 have similar 
minimum line widths (Av = 0.21 kms -1 in Bl and Av = 0.30 kms -1 in B2). Line widths in Oph 
B3 are extremely narrow, with a mean Av = 0.27 kms -1 , a minimum Av = 0.21 kms -1 and an 
rms variation of only 0.04 kms -1 . 

Two locations in B2, both near Class I protostars, show intriguing gradients in line-of-sight 
velocity and line width, which we show more closely in Figure [TUJ Figure [TOa shows the first, 
at a location coincident with the millimeter continuum object B2-MM8, which is found at the 
thermal dus t emission peak in B2 a nd is the most massive clump identified in Oph B (M ~ 1 Mq , 
Motte et al.1 : 1 Johnstone et al I booobh . As noted previously, MM8 is associated with neither an N2H + 



clump nor a maximum of integrated N2H + 1-0 intensity. The nearby Class I protostar, Elias 32, 
may be associated with MM8, but is also offset from the clump peak by ~ 15" to the southeast. 
Across the ~ 30" width of the clump (~ 3600 AU at 120 pc), v hSR varies from ~ 3.8 kms -1 in the 
west to ~ 4.2 kms - 1 in the east. We find n o local decrease in line width as might be expected 



for a starless clump (IDi Francesco et al.ll2004l ). Instead, Figure [TUb shows how at this location Av 
increases from east to west, from Av ~ 0.5 kms -1 to Av ~ 0.9 kms -1 . The resulting velocity 
gradient, given the ~ 3600 AU diameter of the clump, is ~ 23 kms -1 pc -1 . The variation in line 
width in particular appears centered on the nearby protostar. Alternatively, we noted in Paper I 
a region of wide NH3 line width located between Oph B2 and B3, which we suggested was due to 
the overlap of the two Cores with different v LSR along the line of sight. In this case, the wide line 
widths seen near MM8 also may be due to this overlap. 

Figures [TUb and [TQH plot v LSR and Av at the second intriguing gradient, near the Class I pro- 



tosta r Elias 33 in Oph B2. Elias 33 is associated w ith the millimeter object B2-MM10 (jMotte et al 



19981 . also 162729-24274 in Ijprgensen et al.ll2008l ). Again, no N2H + clump was found coincident 



with the protostar. Both v LSR and Av show strong variations at the protostellar location. An 
abrupt change in v LSR from ~ 4.1 kms -1 to ~ 3.7 kms -1 is seen in over a single 18" beam width 
(~ 2200 AU) from west to east. The resulting velocity gradient is ~ 38 kms -1 pc -1 . A narrow 
ridge of large line width, Av > 1.0 kms -1 , is found with an elongation perpendicular to the t> LSR 
gradient. On either side, line widths are narrow (Av < 0.6 kms -1 ). 
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4-1.2. A^i? Density, Column Density and Fractional Abundance 



Given the N2H + 1-0 line characteristics Av, Tt t and T ex determined through HFS fitting, and 
again assuming LTE, we can calculate the volume density, n ex , and the N2H + column density, 
7V(N 2 H+), in Oph B. 



The volume density is calculated following ICaselli et al.l (|2002al ) using two-level statistical 
equilibrium: 



flex 



bg 



(hv/k)(l-T ex /T k ) 

where the critical density, corrected for trapping, of the N2H + 1-0 line is given by: 



(1) 



n„ 



1 - exp(-r h 



(2) 



' hf 



The opacity T h} is the opacity of a typical hyperfine component. For N2H + 1-0, we take r hf = r/9. 
We assume that the N2H + kinetic temperature is equal to the Tk derived in Paper I from NH3 
observations. This assumption introduces some uncertainty in the resulting n ex on the order of 
~ 10 — 50 % if a difference in Tk exists of ~ 2 — 5 K between gas traced by NH3 and N2H + emission 
(a temperature difference > 5K is unlikely given the mean Tk = 15 K for Oph B determined in 
Paper I and the factor ~ 10 — 100 difference in the expected gas densities traced by NH3 and 
N2H + ). Mean, standard deviation, minimum and maximum densities for Oph Bl, B2 and B3 
are given in Table For Oph B as a whole, (n ex ) = 2.4 x 10 5 cm -3 , with an rms variation of 
2.5 x 10 5 cm -3 . Separately, the mean n ex for Oph Bl and B2 (2.0 x 10 5 cm -3 and 3.1 x 10 5 cm -3 , 
respectively) are lo wer by factors of ~ 2 and ~ 3.5 than the mean volume densities calculated by 
Motte et al.l (|1998l ) from 1.3 mm continuum observations. We find only moderate variation in n ex 
across the Core, and consequently the densities determined for the N2H + clumps are similar to the 
mean Core values. This result is discrepant wi th the n(H2) ~ 10 7 cm -3 densities for the small-scale 



continuum objects in Oph B, determined bv lMotte et all but agrees within a factor of ~ 2 with 
the average 850 um c l ump densities (which we calculated from reported clump masses and radii) 
from I Jorgensen et al.1 (|2008l ). 



We follow |Pi Francesco et al.1 (|2004l ) and derive the total N2H + column density A r (N2H + ) 



using 



3/i 1 0r 



tot 



8tt 3 h 2 2y/\^2 {J + 1) exp(-Ej/kT ex ) \ hB 3 J 1 - exp(-hv/kT ( 



J\iTp' 



+ 



Avt 



(3) 



where the l ower rotational ley el number J = 0, the dipole moment of the N2H + molecular ion 
He = 3.4 D (jGreen et al.lll974l ). the N 2 H+ rotational constant B = 46.586702 GHz (jCaselli et al. 
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19951 ) and the rotational level energy above ground of the lower level Ej = J(J + l)hB. Here, 
the rotational partition function Q ro t has been approximated by its integral form for a linear 
molecular ion, Q rr >t = kT ex /hB + 1/3. Note that this calculation is equivalent to that presented by 



Caselli et all (|2002bl ). Figure Et shows the A(N 2 H+) distribution in Oph B. 



We also calculate the H2 column density, N(H.2), per pixel in Oph B from 850 /im continuum 
data using 



A(H 2 ) = S u /[n m fim H K u B v (T d )], 



(4) 



where S v is the 850 fim flux density, Q m is the main-beam solid angle, /i = 2.33 is the mean molecu- 
lar weight, m H is the mass of hydrogen, k u is the dust opacity per unit mass at 850 fim, and B v (Td) 
is the Pl anck function at the dust temperature, T^. As in P aper I, we take k v = . 018 cm 2 g , 
following IShirlev et al.l ((2000), using the OH5 dust model from lOssenkopf fc Hennina ()1994l ) which 
describes grains that have coagulated for 10 5 years at a density of 10 6 cm -3 with accreted ice man- 
tles, and incorporates a gas-to-dust mass ratio of 100. The 15" resolution continuum data were 
convolved to a final angular resolution of 18" and regridded to match the Nobeyama N2H" 1 " data. 

We assume that the dust temperature per pixel is equal to the kinetic temperature re- 
determined in Paper I for NH3 in Oph B, where we found (Tk) = 15 K. This assumption is 
expected to be good at the densities probed by the NH3 data (n > 10 4 cm~ 3 ), when thermal 



coupling between the gas and dust by collisions is expected to begin (jGoldsmith fc Langerlll978l ). 
Our analysis in Paper I, however, suggested that NH3 observations may not trace the densest 
gas in Oph B. As we expect the denser gas to be colder in the absence of a heating source, this 
assumption may systematically overestimate in the Core. The derived N(H.2) values would then 
be systematically low. There is additionally a ~ 2 0% uncertainty in the continuum flux values, and 
estimates of k„ can vary by ~ 3 (jShirley et al.ll2000l ) . Our derived H2 column densities consequently 
have uncertainties of factors of a few. 

Due to the chopping technique used in the submillimeter observations to remove the bright 
submillimeter sky, any large scale cloud emission is necessarily removed. As a result, the image 
rec onstruction technique p roduces negative features around strong emission sources, such as Oph 
B (1 Johnstone et al.l l2000al ). Emission on scales greater than 130" is suppressed by convolving 
the original continuum map with a a = 130" continuum beam, and subtracting the convolved 
map from the original data. While the flux density measurements of bright sources are likely 
accurate, the column densities measured at core edges may be underestimates. Following Paper I, 
we thus limit our analysis to pixels where S u > 0.1 Jybeam -1 , though the rms noise level of the 
continuum map is ~ 0.03 Jybeam -1 . For a dust temperature = 15 K, this flux level corresponds 
to A(H 2 ) ~6x 10 21 cm -2 . 

Using A(H2) and A(N2H + ), we next calculate per pixel the fractional abundance of N2H + 
relative to H2, A(N2H + ), in Oph B. In Figure [9H, we show the distribution of A(N2H + ) in Oph B. 
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The mean, rms, minimum and maximum values of the derived N2H + column density and fractional 
abundance are given in Table [5j while specific values for individual N2H + clumps are given in Table 
HJ The iV(H2) and consequently the X(N2H + ) uncertainties given in Table H] include the ~ 20% 
uncertainty in the submillimeter continuum flux values only; uncertainties in and k v are not 
taken into account. 

The mean N2H + column density (iV(N2H + )) = 5.5 x 10 12 cm~ 2 in Oph B. As seen in Table O 
N2H + column densities in Oph B2 are slightly greater, on average, than in Oph Bl. Bl, however, 
contains the largest N2H + column density peak, with a maximum iV(N2H + ) = 1.7 x 10 13 cm~ 2 
found near (but not at the peak position of) the N2H + clump B1-N4. The smallest average column 
density is found towards B3, and the maximum value in B3 is ~ 2 — 4 times less than peak values 
found in both Bl and B2 (iV(N2H + ) = 1.1 x 10 13 cm -2 towards B2-N9; similar high values are also 
found towards B2-N2). With the exception of the Bl iV(N2H + ) peak, observed peaks in A r (N2H + ) 
are only factors of ~ 2 above the mean values. 

In Figure [9b, we show the locations of N2H + clumps identified in ^3.1.11 over the A r (N2H + ) 
distribution in Oph B. The peak emission locations of N2H + clumps are often found near, but are 
not perfectly coincident with, local 7V(N2H + ) maxima. Clumps identified in continuum emission 
have varying correspondence with A r (N2H + ) maxima. In Oph Bl, two of three 850 /im continuum 
clumps are coincident with the two largest iV(N2H + ) maxima, while in Oph B2 only two of five 
continuum objects is co-located with an iV(N2H + ) maximum. The embedded protostars in B2 are 
not associated with A r (N2H + ) maxima. 

The mean fractional N2H + abundance for Oph B, (X(N2H + )) = 8.0 x 10~ 10 , is greater by 
a factor of ~ 2 — 4 than previous measurements in isolated starless cores 



Tafalla et al. 


2002, 


2006; 


Keto et al. 


2004) 



Bl is greater by a factor ~ 2 than in Oph B2, although the maximum values are similar. The 
greatest N2H + abundances are found towards the sub-Core edges and the connecting filament. 
Conversely, the smallest X(N2H + ) values are located at the centres of Bl and B2. Little small- 
scale variation is found in the sub-Core interiors. We discuss further in §5.2 the general trend in 
X(N2H + ) with H2 column density. 

At the peak locations of N2H+ clumps in B2, the mean X(N2H + ) = 3.2 x 10~ 10 is slightly 
greater than typical values found in isolated starless cores. N2H + clumps, however, are not located 
at clear X(N2H + ) maxima at the 18" resolution of the single-dish data, as shown in Figure [9ji 
where we plot the N2H + clump locations over the X(N 2 H+) distribution in Oph B. In Oph B2, 
two of five continuum clumps are found towards N2H + abundance minima, and no clear abundance 
minima or maxima are associated with the remaining three. In Oph Bl, two of three continuum 
clumps are co- located with X(N2H + ) minima, while a good abundance measurement was not found 
at the position of the third. 



-16- 



4.2. Interferometer Results 

We next look at the line fitting results in the ATCA pointings. Complex line structures were 
frequently found in the interferometer data (see ^3.2|) . and in combination with relatively low S/N 
in many locations a full 7-component Gaussian N2H + 1-0 line fit did not converge. Instead, we 
fit a single Gaussian to the isolated N2H + 1-0 component to determine the line v LSR and Av. In 
Figures 1H[5] and [6] (parts b and c) we show -u LS r and Av over the B2-MM8 region, the eastern N2H + 
object B2-N5, and the south-eastern edge of Oph B2, respectively. While emission was found in 
the north-eastern mosaic (see Figure ED, the relatively low line strengths compared with the mosaic 
rms and complex line structures precluded good fits to the data. 

In general, we find the structures visible in the interferometer data have similar v LSR within 
uncertainties, and universally narrower Av by a factor ~ 2, compared with the single-dish data 
at the same location. This difference can be seen qualitatively in Figure [Til where we plot the 
Nobeyama and ATCA spectra towards B2-MM8, B2-N5, B2-N6, B2-N7 and B2-N8. (We show 
the B2-N6, B2-N7 and B2-N8 spectra at the interferometer data peak, as the ATCA data show 
integrated intensity maxima offset by ~ 5" from the single-dish data towards B2-N6 and B2- 
N8, while significant interferometer emission was not seen towards B2-N7.) The B2-N5 v hSR is 
blueshifted by 0.1 kms" 1 in the ATCA data relative to the Nobeyama data, but Figure [TT1 shows 
this is likely due to strong negative features in the interferometer spectrum. 

We find a velocity gradient of ~ 0.2 kms -1 over ~ 10" (0.005 pc, or ~ 33 kms^pc" 1 ) with 
a position angle of ~ 75° across the small interferometer peak north of B2-N6 (see Figure Eb). If 
this gradient is associated with rotation of the clump, this value is large compared with typical 
gradient values observed at low resolution, but is compara ble to high resolutio n N2H + 1-0 results for 



dense starless and protostellar clumps. In Ophiuchus A, [Di Francesco et all (J2004J) found velocity 



gradients of 0.5 kms 1 o ver 0.01 pc, or 50 kms 1 pc , towards several submillimeter clumps. In 



addition. iLee et all (|2007l ) found a velocity gradient of 30 / cos(i) kms" pc~ ) within 30" (0.04 pc 
at 300 pc) towards the L1521B protostellar group, perpendicular to an observed outflow direction 
and suggestive of fast rotation. The gradient seen in L1521B is both an order of magnitude larger 
than that seen at lower resolution and in the opposite direction (a neighbouring core was in cluded 



in the low resolution data). In the protostellar core IRAM 04191, iBelloche Sz Andrei (J2004J) found 
a mean velocity gradient of ~ 17 kms -1 pc -1 through interferometer N2H + observations, which is 
a factor ~ 5 larger than that seen in single-dish data. 

B2-N6 is also associated in the interferometer data with extremely narrow N2H + 1-0 line 
widths, with a mean Av = 0.21 ± 0.04 kms" 1 averaged over an 8" x 5" area matching the beam 
FWHM. The observed single-dish line width at the same location is 0.39 ± 0.01 kms" 1 , listed in 
Table HI These are the narrowest lines found in the Oph B2 Core. (Slightly greater line widths are 
found in Oph Bl towards the B1-N3 and B1-N4 clumps when fitted by a two-velocity component 
structure in the single-dish data, see $H) The velocity gradient, small size and narrow line widths 
associated with B2-N6 in the interferometer data make it an interesting emission feature in B2, 
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which we discuss further in §5.2.21 



5. Discussion 

5.1. General trends 

We show in Figure [12] the distribution of the ratio of the non-thermal dispersion to the sound 
speed, o" NT /c s , n ex , 7V(N2H + ) and X(N2H + ) with iV(H2) in Oph B, omitting pixels where the 
850 fim continuum flux S v < O.lJybeam -1 . Data points represent values for 18" pixels, i.e., 
approximately the Nobeyama beam FWHM. We show individually pixels in Oph Bl and B2, 
and also show peak va lues for N2H + clumps, identified in §3.1.11 and 850 fim continuum clumps 



(jJorgensen et al.l 120081 ) . Note that no Oph B3 values are plotted due to the lack of continuum 
emission at the B3 location. The Figure shows that the N2H + clumps do not reside at the highest 
H2 column densities, but scatter over the range of -/V(H2) calculated above our submillimeter flux 
threshold. The submillimeter clumps tend to be found at higher 7V(H2) values than the N2H + 
clumps, on average, but also show a spread in peak iV(H2). 

In Paper I, we used the ratio of NH3 (1,1) and (2,2) emission lines in Oph B to determine the 
kinetic temperature Tk in each pixel, which we can use to calculate <r NT and c s for N2H + emission 
in Oph B. Given Tk, Cnt = yj o% s — ^bTk / (fJ-mol^H) > where ks is the Boltzmann constant, m H is 
the mass of the hydrogen atom, fi mo i is the molecular weight in atomic units (/% 2H + = 29.02) and 
dobs = Av I (2\/21n2). We list values and propagated uncertainties for individual N2H" 1 " clumps in 
Table HI and give the mean, rms, minimum and maximum c NT / c s in Table [5] for each of Oph Bl, 
B2 and B3. 

We find a mean <r NT / c s = 1.02 for N2H + 1-0 emission across Oph B with an rms variation of 
0.49, indicating that the non-thermal motions are approximately equal, on average, to the sound 
speed. When looked at separately, we find a difference in the <r NT / c s ratio between Oph Bl and 
B2. In Oph Bl, the mean a NT / c s = 0.71, with an rms variation of 0.16. In Oph B2, the mean 
°nt / c s = 1-26 with an rms variation of 0.4. Oph B3 is dominated by nearly thermal line widths, 
with the mean a NT / c s = 0.39 and an rms variation of only 0.13. 

We find little variation in cr NT / c s with iV(H2) in Oph B, as shown in Figure [T2"a, The difference 
in the mean a NT / c s ratio between Bl and B2 is clear, as is the significantly larger scatter of non- 
thermal line widths in B2 than in Bl. The scatter in cr NT / c s extends to the highest H2 column 
densities. In Bl, all non-thermal line widths a NT ^ c s , shown by the dashed line. Similar results are 
found for most N2H + clumps, and approximately half the continuum clumps. At 18" resolution, 
the observed non-thermal line widths do not extend to arbitrarily low values, as the data show a 
cutoff in the a NT / c s ratio at cr NT / c s ~ 0.5 for both Bl and B2, below which no data points are 
found (but note that in §4.21 we found narrow line widths at high resolution). For Tk = 15 K, the 
mean temperature found in Oph B, the cj^t 

/ c s cutoff corresponds to a NT ~ 0.12 km s 1 . Note 
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that the velocity resolution of the N2H + data, Av r 
observed lower a NT / c s limit. 



0.05 kms , is significantly less than the 



In Paper I, we discussed the possible origins of the large non-thermal motions in Oph B, 
and based on timescale arguments ruled out primordial motions (i.e., motions inherited from the 
parent cloud). It is interesting that most of the gas traced by N2H + 1-0 in Oph B2 remains 
dominated by transonic non-thermal motions, while the line widths in Oph Bl reveal the Core is 
significantly more quiescent at the ~ 10 5 cm" 3 densities traced by N2H + 1-0 emission. The most 
obvious so urce of this differ e nce is the presence of embedded protostars in Oph B2 while Oph Bl 
is starless. iKamazaki et al.l (|2003l ) found an outflow in CO 3-2 emission in Oph B2 centred near 
Elias 33 and Elias 32 and oriented approximately east- west, but were unable to pinpoint which 
protostar was the driving sour ce. In upcoming results from the JCMT Gould Belt Legacy survey 
(|Ward-Thompson et al.ll2007bh . which mapped Oph B in CO 3-2, 13 CO 3-2 and C 18 3-2 at 14" 
resolution, this highly clumped outflow is shown to extend over more than 10' (i.e., 0.4 pc at 120 pc, 
White et al. 2009, in preparation). The outflow axis is also aligned such that the outflow does not 
appear to be impacting Oph Bl significantly, potentially explaining the difference in non-thermal 
line widths between the two sub-Cores^|. 

N2H + is not generally thought to be a tracer of protostellar outflows, since it is expected to 
be destroyed quickly through reactions with CO, which evaporates from dust grains in the higher 
temperature gas near the driving protostar. On small scales in B2, however, some variations in v LSR 
and Av appear correlat ed with the presence of nearby Class I protostars (see §4.1.11 and Figure 

Co 



Chen et al.1 (|2008l ) suggest N2H + can be entrained in protostellar jets before a molecular 
outflow releases CO from grain sur faces which then can des troy N2H + . In a small N2H + 1-0 



survey of the Serpens NW cluster, IWilliams k, Mversl (|2000l ) found that the largest <r NT values 
(<t nt > 0.6 kms -1 , greater than seen here in Oph B2) occurred in cores containing protostellar 
sources which power strong outflows. If caused by the substantial protostellar outflow, the large 
line widths on a global scale in Oph B2 suggest that protostellar outflows are able to inject additional 
turbulence into the high density gas in cluster- forming Cores, thereby increasing the mass at which 
clumps become gravitationally unstable and altering the consequent fragmentation and evolution 
of existing embedded clumps. iDi Francesco et al.1 (120011 ) did not see this effect in the N2H + 1-0 
observations of the circumstellar dense gas associated with the protostellar, outflow-driving source 
NGC 1333 IRAS4B, suggesting that the ability of the outflow to inject turbulent energy may be 
determined by additional factors, such as the collimation or strength of the outflow. A comparison 
of the outflow properties in these regions would be useful to probe further the impact of outflows 
on cluster forming dense gas. 

Alternatively, the larger line widths in B2 may be due to global infall motions. In Paper I, 
we determined a virial mass M v i r ~ 8M Q for Oph B2 based on NH3 line widths, which is a factor 
~ 5 less than Core mass estimates based on thermal dust continuum emission (with a factor ~ 2 



2 JCMT Spring 2009 newsletter, |http://www .jach.hawaii.edu/J CMT/publications/newsletter/n30/jcmt-n30.pdf 
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uncertainty) . In contrast, M/M v i r ~ 1 in Oph B l. Some evidence for infall has been o bserved in 



self-absorbed molecular line tracers in Oph B2 ([Andre et al J 120071 : iGurnev et al.ll2008l ). but it is 



difficult to determine unambiguously given the complex outflow motions also present. Additionally, 
the consistent line widths found in NH3 gas for both Bl and B2 (discussed further in §5.3.11 in 
comparison with the N2H + results) suggest a common source for non-thermal motions at gas 
densities n ~ 10 3-4 cm -3 , which then impacts differently the gas at higher densities in the two 
sub-Cores. 

Figure 112b shows the volume density Tl ex cLS Si function of A^B^). Nearly all pixels have 
n ex > 10 5 cm" 3 , and we find slightly greater n ox at higher column densities. The N2H + and 
continuum clumps scatter over the full range of density values found in Oph B. 

We find a slight trend of increasing N2H + column density, A(N2H + ), with iV(H2), shown in 
Figure [TZb. although the mean A(N2H + ) changes by less than a factor of ~ 2 over the order of 
magnitude range in N(H.2) found in Oph B. Several pixels in Oph Bl, and one N2H + clump, have 
significantly greater A(N2H + ) for their N(H.2) values compared with the rest of the data. The 
small gradient in A"(N2H + ) with (H2) leads to a trend of decreasing fractional N2H + abundance 
with increasing H2 column density, shown in Figure fl~2"H . We find N2H + abundances decrease by 
an order of magnitude between the limiting N(H.2) threshold and the peak N(H.2) values. We fit a 
linear relationship to logX(N2H + ) versus log A^H^) for Oph Bl and B2 separately. In Bl alone, 
where we have relatively few data points, we find log A(N2H + ) is consistent within the uncertainties 
with being constant with A^Efe), but in conjunction with the Oph B2 data, we find 

log A(N 2 H+) = (7.1 ± 0.9) - (0.74 ± 0.04) x log A(H 2 ) (5) 
using the same H2 column density limits as in Figure H21 



5.2. Small- Scale Features 

We next discuss the physical properties of the small-scale structure present in Oph B, including 
the N2H + clumps, continuum clumps, and embedded protostars, and describe in detail the compact, 
thermal object B2-N6. 



5.2.1. Comparison of N2H + clumps, continuum clumps and protostars 



We list in Table [6] the physical parameters d erived from N2H + 1-0 emission to wards 850 um 



conti nuum clump locations (IJorgensen et al.ll2008l ) and embedded Class I protostars ([Enoch et al 



2009). Columns are the same as in Table [H In Table El we list the mean values of each physical 



parameter towards N2H + clumps, continuum clumps and protostars. Note that we were only able 
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to solve for all parameters listed towards a single protostar, and discuss below only those mean 
parameters where we obtained values from three or more objects. 

We find no difference in the mean vlsr between N2H + clumps, continuum clumps and proto- 
stars. There is a clear reduction in Av towards N2H + clumps relative to continuum clumps (mean 
Av = 0.49 kms _1 and 0.68 kms -1 , respectively), but excluding one continuum clump (162712- 
24290) reduces the mean Av to 0.58 kms -1 . Protostars are associated with wider mean Av, but 
the larger mean is also driven by large Av towards a single object. The differences between N2H + 
and continuum clumps are marginally significant given the variance of Av across the entire Core is 
only 0.08 kms -1 in Oph Bl and 0.21 kms -1 in Oph B2. We find similar N2H + line opacities and 
excitation temperatures towards continuum and N2H + clumps. Continuum clumps are associated 
with larger non-thermal motions, with a mean <jnt/cs = 1-24 which is 1.5 times that associated 
with N2H + clumps (mean ctnt/cs = 0.86). N2H + column densities are similar towards N2H + 
and continuum clumps (within ~ 25 %), but N2H + clumps are associated with smaller iV(H2) and 
hence greater X(N2H + ) by a factor > 2. In fact, the mean N2H + abundance towards continuum 
clumps is less than X(N2H + ) averaged over the entire Oph B Core. Derived volume densities are 
the same for N2H + and continuum clumps. 

In summary, the N2H + clumps are associated with smaller line widths and subsonic non- 
thermal motions relative to continuum clumps and protostars, but the difference in mean values 
is similar to the variance in line widths across the Core. The most significant difference between 
N2H + and continuum clumps are found in their N2H + abundances, with a mean X(N 2 H+) towards 
continuum clumps that is lower than AT(N2H + ) in N2H + clumps. 



5.2.2. Oph B2-N6 



The B2-N6 clump is a pocket of narrow N2H + 1-0 line emission within the larger, more 
turbulent B2 Core, with a large fraction (> 50%) of the emission at the clump peak coming from 
the small size scales probed by the ATCA observations (see Figure [5]). Although we do not consider 
all N2H + clumps to trace underlying physical structure in Oph B, the remarkable properties of 
B2-N6 suggest that further attention be paid to that particular clump. Line widths measured 
from the interferometer data are a factor ~ 2 less than those recorded from the single-dish data 
(Av = 0.21 ± 0.04 kms -1 compared with 0.39 ± 0.01 km s" 1 ). The millime ter continuum object 



B2-MM15, identified through a spatial filtering technique ( Motte et al. 19981). lies offset by ~ 12" 



from the peak interferometer integrated intensity (see Figure [5]) . iMotte et al.l calculate a mass of 
only 0.17 M associated with the B2-MM15 clump, but based on its small size (unresolved at 15" 
resolution, using an estimate of r ~ 800 AU) derive an average density of n ~ 1.2 x 10 8 cm -3 . This 
density is significantly larger than the densities calculated in §4.1.21 at 18" resolution, but higher 
sensitivity N2H + 1-0 interferometer observations are needed to probe the clump density on small 
scales. 
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We calculated in Paper I a gas kinetic temperature Tk = 13.9 ± 1.1 K for the co-located 
B2-A7 NH 3 (1,1) clump. We find cr NT = 0.06 kms' 1 for B2-N6, equal to the expected N 2 H+ 
thermal dispersion, and a NT /c s = 0.26. A simi lar compact N2H" 1 " 1-0 c lump with nearly thermal 
line widths was found in Oph A (Oph A-N6, |Pi Francesco et al.1 12004 ) . While both clumps are 
associated with larger-scale N 2 H + emission, these objects are only distinguishable as compact, 
nearly purely thermal clumps when observed with interferometric resolution and spatial filtering. 



Di Francesco et al.l proposed Oph A-N6 may be a candidate thermally dominated, critical 
Bonnor-Ebert sphere embedded within the more turbulent Oph A Core. Such objects have sizes 
comparable to the cutoff wavelength for MHD waves a nd are in equilibrium between their self- 
gravity and internal and external pressures. iMversI (|1998l ) called these objects 'kernels', and showed 
that they could exist in cluster-forming cores with FWHM line widths Av > 0.9 kms -1 and column 
densities A^(H 2 ) > 10 22 cm -3 , with sizes comparable to the spacing of protostars in embedded 
clusters, suggesting that a population of kernels within a dense Core could form a stellar cluster. 
The physical conditions in both Oph A and Oph B2 fit these requirements. The B2-N6 clump, 
however, is significantly smaller (0.004 pc compared with a predicted ~ 0.0 3 pc) an d less massive 
(O.17M compared with a predicted ~ 1M ) than the objects discussed by iMvers . 



If we calculate the mass of a critical Bonnor-Ebert sphere given the ATCA N2H + 1-0 line 
width and a radius of ~ 800 AU, where M/R = 2Aa 2 /G feonnoilllflSfih. w e find M BE = 0.02 M 



0; 

Motte et al. We find a virial mass 



or approximately 10 x smaller than the mass determined by 
M v i r = 0.04 Mq, assuming B2-N6 is a uniform sphere and M, 
i s acc urate, this suggests that while small, B2-N6 may be gravitationally unstable. lAndre et al 



5a 2 R/G. If the clump mass 



(|2007l ) made a tentative detection of infall motions towards the clump through observations of 
optically thick lines such as CS, H2CO or HCO + . Thus B2-N6 may be at a very early stage of 
clump formation, and may eventually form an additional low mass protostar to the three YSOs 
already associated with dense gas in Oph B2. 



5.3. Comparison of N2H + and NH3 emission in Oph B 

We next compare the physical properties of the gas in Oph B derived here from N2H + observa- 
tions with those derived from NH3 observations, described in Paper I. In the following discussion, 
all comparisons are made after convolving the NH3 data to a final angular resolution of 18" to 
match the Nobeyama N2H + observations. 



5.3.1. v LSR and a NT / c s 

We first look at the velocities and non-thermal line widths of the gas traced by N2H + gas 
and those determined from NH3 emission in Paper I. Figure [T3"a shows the distribution of v LSR 
derived from N2H + and NH3 emission in Oph B. There is no significant difference in the measured 
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line-of-sight velocity between the two dense gas tracers, with the mean u LSR = 3.98 kms 1 for NH3 
and 4.01 kms" 1 for N 2 H+. 

In Figure 113b . we show the distribution of c NT / c s as determined with N2H + or NH3 using 
the NH^-derived Tk values. The non-thermal line widths measured in N2H + 1-0 are substantially 
smaller than those measured in NH3 in Paper I, and additionally the variation in the relative 
magnitude of the non-thermal motions between Oph Bl and B2 was not found in NH3 emission. The 
mean NH3 <r NT / c s = 1.64 in Oph B, and does not vary significantly between Bl ((<r NT / c s ) = 1.62) 
and B2 ((cr NT /c s ) = 1.68). 

The assumption that Tk is similar for both N2H + and NH3 is reasonable if both molecules 
are excited in the same material. The good agreements between the locations of NH3 and N2H + 
clumps, illustrated in Figure El and the u L sr) illustrated in Figure [T3h . support this assumption. 
The substantial offset in the cr NT / c s ratio seen in Figure 113b . however, between NH3 and N2H + 
emission indicates significantly different motions are present in the gas traced by NH3 than in the 
gas traced by N2H + . The gas densities traced by N2H + are ~ 1 — 2 orders of magnitude greater 
than those traced by NH3 (~ 10 5 cm -3 and ~ 10 3-4 cm~ 3 , respectively). In fact, if we attempt to 
derive the gas kinetic temperature by assuming equal non-thermal motions for the N2H + 1-0 and 
NH3 (1,1) emission (effectively assuming N2H + and NH3 trace the same material), the resulting 
average value is a highly unlikely Tk = 190 K for Oph B2. Given the narrow N2H + lines observed 
in Oph Bl, no physical solution for Tk can be found under the assumption of equal <r NT . 

Starless cores are typically found to be well describe d by gas temperatures th at are either 



constant or decreasing as a function of increasing density (|Di Francesco et al.l 120071 ). Given that 
N2H + traces denser gas, we then expect Tk from NH3 measurements to be biased high in starless 
cores, e.g., gas traced by N2H + should be colder than gas traced by NH3. The mean Tk found in 
Oph B in Paper I was 15 K. With a lower Tk, the sound speed would be smaller and the returned 
cr NT would be larger on average. A gas temperature Tk = 10 K rather than 15 K would increase 
the average N2H + <r NT / c s ratio in Oph B by ~ 20 — 25%. This increase, while significant, is 
not large enough for the mean N2H + <r NT / c s to match the NH3 results in Oph B. With both 
a constant temperature or decreasing temperature with density, non-thermal motions in Oph Bl 
and B3 would remain subsonic, while Oph B2 would still be dominated by transonic non-thermal 
motions. Alternatively, it is possible that the two Class I protostars in Oph B2 could raise the 
temperature of the dense gas above that traced by NH3 (no Tk difference was found between 
protostars and starless areas in Paper I). In this case, the returned N2H + <r NT would be smaller, 
and would further increase the differences in magnitude of non-thermal motions seen between the 
NH3 and N2H + emission in Oph Bl. In B2, a temperature of 20 K would decrease the mean 



g NT / c. q ratio to 1.1. In a study of dust temperatures in the Ophiuchus Cores, IStamatellos et al 



(|2007l ) showed, however, that embedded protostars in the Cores will only heat very nearby gas and 
will not raise the mean temperature of the Cores by more than ~ 1 — 2 K, so an average Tk = 20 K 
over all Oph B2 is unlikely. 
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5.3.2. Double peaked line profiles 

In Figure [T-41 we show the isolated F\F — ► F[F' = 01 — ► 12 component of the N2H + 1-0 
emission line towards the three N2H + clumps in Oph Bl (Bl-Nl, B1-N3, and B1-N4) which show 
double peaked line profiles, as described in SJH Figure HH also shows a single Gaussian line profile 
overlaid on each N2H + spectrum, which represents the v LSR and Av of the 18 component Gaussian 
fit to the NH3 (1,1) emission at the clump peak, normalized to a peak amplitude of 1 K. We plot the 
Gaussian fit for clarity since there are no isolated components in the NH3 (1,1) hyperfine structure. 
Additionally, spectra of C2S 2i — 1q emission line at the peak clump locations are shown. Note 
that over all Oph B, C2S emission was only detected towards the southern tip of Oph Bl, and 
thus no significant C2S emission was found at the peak location of Bl-Nl. The C2S emission was 
observed with the GBT at ~ 32" spatial resolution and 0.08 kms -1 spectral resolution, with the 
observations and analysis described in Paper I. Arrows show the locations of the fitted v hSR for all 
species, including both N2H + velocity components. 

The clumps which show the double-peaked line profile structure are the clumps with the 
highest total N2H + 1-0 line opacities based on a single velocity component fit, with r = 5,6 and 
9, respectively for Bl-Nl, B1-N3 and B1-N4. The F\F -> F[F' = 01 ^ 12 component is expected 
to have an opacity Tj so = l/9r, so Ti so = 0.6,0.7, and ~ 1 for the three clumps, based on a 
single component fit. If the line is optically thick, however, this fit is likely to underestimate the 
true opacity given the missing, self-absorbed flux. This is suggestive that the lines are indeed 
self-absorbed. Since the NH3 and C2S emission is optically thin (Paper I), however, if the N2H + 
emission was self-absorbed due to high optical depth we would expect the emission peak of both 
NH3 and C2S to be found between the two N2H + components in v LSR , which is not the case. 
In all three clumps, we find the t> LSR of the NH3 emission more closely matches the v LSR of the 
red N2H + component, while the C2S emission in two clumps (v LS r = 3.67 kms -1 ) more closely 
matches the v LS n of the blue N2H + component (v hS n = 3.85 kms -1 and 3.62 kms -1 for B1-N3 and 
B1-N4, respectively). This behaviour suggests the N2H + 1-0 emission is not sel f-absorbed t o ward s 



these positions. A similar offset of C2S emission relative to N2H + was found bv lSwift et al.l ((2006) 
towards the starless Core L1551, however the case of L1551 the C2S emission is redshifted with 
respect to the systemic velocity of the Core rather than blueshifted, as we find in Bl. 

The GBT NH3 (1,1) observations described in Paper I found extensive emission around Oph 
B, such that it is possible to determine the v LSR and Av of NH3 beyond where continuum emission 
traced the Oph B Core. We use an intensity threshold of 2K in the NH3 (1,1) main component 
to separate Core and off-Core gas. The off-Core gas surrounding Oph Bl has a mean v LSn = 
3.72 kms -1 with an rms variation of 0.13 kms -1 , which is significantly different from the mean 
Bl Core v hSR = 3.96 kms -1 and rms variation of 0.15 kms -1 . Both the C2S emission and the 
blue N2H + component are thus more kinematically s imilar to the off-Core NH3 gas. The critica l 



densities of C2S 2i — lo and NH3 (1,1) are similar (jSuzuki et al.l 119921 ; iRosolowskv et al.ll2008l ). 
Kinematically, we find <r NT / c s = 0.7 for the C2S emission if we assume the Tk determined from 
NH3 emission accurately represents the temperature of the gas traced by C2S. This result closely 
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matches the N2H + results in Oph Bl, but the non-thermal motions in the gas traced by NH3 are 
significantly larger (£ )5.1I and Paper I). It is impossible to match the non-thermal C2S and NH3 
motions for a given Tk, since the NH3 (1,1) line cr NT > a b s of the C2S. 

At densities n > a few xlO 3 — 10 4 cm" 3 , chemical model s predict C?S is quickly depleted from 



the gas phase (timescale td ep ~ 10 5 yr; Millar Sz Herbst 199C ; Bergin &: Langer 19971), and accord- 
ingly the molecule is obse rved to be a sensitive tracer of depletion in isolated cores (|Lai &: Crutcher 
2000; iTafalla et alj 12006). The volume density in Oph Bl is greater than that required for C2S 



depletion (n > 10 5 — 10 6 cm -3 from £ 14.1.21 and continuum observations), so we would not expect 
to see any C2S emission unless the gas has only been at high density for t < t^ v . C2S was only 
detected towards southern Bl. An explanation, therefore, for both the presence of C2S in southern 
Bl and for its velocity offset relative to the dense gas is that gas from the ambient molecular cloud 
is accreting onto Bl, reaching a density high enough to excite the C2S emission line. This material 
would be chemically 'younger' than the rest of the Bl Core gas, such that the C2S has not had 
enough time to deplete from the gas phase. A similar result was fo und in a multi-species study 
of the starless L1498 and L1517B cores, where ITafalla et ajj (J2004J) conclude that non-spherical 
contraction of the cores produced asymmetric distributions of CS and CO, where the CS and CO 
'hot spots' reveal the distribution of recently accreted, less chemically processed dense gas. 



A single 850 /im continuum clump (162715-24303; iJorgensen et al.l l2008l ) was identified in 
southern Oph Bl where we find B1-N3 and B1-N4. The authors found a clump mass M = 0.2 Mq, 
which is a factor ~ 2 less than the virial mass M V i r = 0.4M Q we calculate based on the mean 
line width (of the single velocity component fit) of B1-N3 and B1-N4 (Av = 0.52 kms , or 
a = 0.22 kms -1 ) and the continuum clump radius. This suggests that the clump is currently stable 
against collapse. If the clump is gaining mass through ongoing accretion from the ambient gas, 
however, then eventually the non-thermal support may not be large enough to prevent gravitational 
collapse, leading to the formation of a low-mass protostar. 



5.3.3. Opacity and Excitation Temperature 

We show in Figure [T3b and [13b the distribution of total line opacity r and excitation tem- 
perature T ex derived from hyperfine line fitting of NH3 (1,1) and N2H + 1-0 in Oph B. We find 
higher total opacities in N2H + emission ((r) = 2.5) than in NH3 emission ((r) = 1.0) in Oph B. 
The N2H" 1 " emission is characterized by lower excitation temperatures than those found for NH3 
emission in Paper I, with a mean T ex = 7.1 K for N2H + 1-0 compared with T ex = 9.5 K for NH3 
(1,1). 
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5.3.4- Fractional Abundances and Chemical Evolution 

We next compare the column densities of NH3 and N2H + towards the Oph B Core. We only 
show results for pixels which have well-determined values for both A^NFFj) and iV(N2H + ). In 
Paper I, we found a trend of decreasing fractional NH3 abundance with increasing iV(H2) for both 
Bl and B2. After convolving the NH3 data to match the single-dish N2H + resolution, we find the 
following linear relationships between logX(NH3) and logiV(H2): 



Bl : logX(NH 3 ) = (5 ±4) - (0.6 ± 0.2) log iV(H 2 ) 
B2 : logX(NH 3 ) = (3.5 ± 2.2) - (0.5 ± 0.1) log iV(H 2 ) 



(6) 
(7) 



The slopes found for Oph Bl and B2 are both negative and similar to that found in Equation [5] 
for the N 2 H + fractional abundance as a function of N(H.2). This result is suggestive of a trend of 
increasing depletion of both NH3 and N2H + with increasing N(H.2) in Oph B. 

We show in Figure [15] the NH3 column density, iV(NH3), versus the N2H + column density, 
Y(N 2 H+) in Oph Bl and B2, and also the ratio of A(NH 3 ) to iy(N 2 H+) as a function of A(H 2 ), 
calculated in §4.1.21 (we do not show results in B3 due to the small number of pixels with both 
well-determined NH3 and N2H + column densities, but note mean values below). Note that the 
column density ratio, A(NH3) / iV(N 2 H + ), is equivalent to the ratio of fractional abundances, 
X(NH 3 ) /X(N 2 H+). Figure QBt shows a significant difference in the iY(NH 3 ) / A r (N 2 H + ) ratio 
between Oph Bl and B2. We find the mean and the standard deviation of this ratio are each 
twice as large in Bl (iV(NH 3 ) / 7V(N 2 H+) = 135, a = 52) as in B2 (JV(NHa) / A(N 2 H+) = 65, 
a = 31) and B3 (iV(NH 3 ) / iV(N 2 H+) = 80, a = 56). In Oph Bl and B2, Figure illustrates 
that this difference is largely driven by lower A^(N 2 H + ) values towards Bl, while the spread in the 
iV(NH3) / iV(N2H + ) appears due to a wider spread of A(NH3) values in Bl. We find no variation 
in the A(NH 3 ) / A(N 2 H+) ratio as a function of JV(H 2 ). 



In the clustered star forming region IRAS 20293+3952, iPalau et al.l (|2007l ) found strong 
NH3 and N 2 H + differentiation associated with the presence or lack of embedded YSOs. Low 
iV(NH 3 ) / A r (N 2 H + ) ~ 50 values were found near an embedded YSO cluster, while higher A(NH 3 ) / A(N 2 FF 
300 were found towards cores with no associated YSOs. In a stud y of two dense cores, each contain- 
ing a starless main body and a YSO offset from the core center, iHotzel et al.1 (120041') find a f actor 
~ 2 smaller A(NH3) / X(N 2 H + ) values towards the starless cores compared with lPalau et all with 
abundance ratios (~ 140 — 190, towards the starless gas, and ~ 60 — 90 tow ards the YSO s) similar 
to those found in this study. The relative abundance variation found by IHotzel et al.l is driven 
by a varying NH3 abundance while X(N 2 H + ) remains constant over the cores. Although we find 
both a decrease in N 2 H + column density as well as an increase in NH3 column density drives the 
greater fractional NH3 abundance towards Oph Bl (see Figure PT5|) . both X(N 2 H + ) and A(NH3) 
are greater in Bl, by factors of ~ 2 and ~ 4, respectively, compared with B 2 given the lower H 2 
column densities found in Bl in §4.1.21 In a survey of 60 low mass cloud cores. ICaselli et al.l (|2002al ) 
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find correlations between Af(NH3) and 7V(N2H + ) in starless cores, but the relative column density 
values do not vary significantly between starless and protostellar objects. 

Based on this work and the studies described above, it appears that the relative fractional 
abundance of NH3 to N2H + remains larger towards starless cores than towards protostellar cores 
by a factor of ~ 2 — 6 in both isolated and clustered star forming regions. We find a 1-a positive 
trend in the relative NH3 to N2H + abundance in Oph B2 with increasing distance to an embedded 
protostar, suggestive of a direct impact by the protostars on the relative abundance s. In their 
work on x-ray emission from young protostars in Ophiuchus, I Casanova et al.l (|1995l ) show that 
protostellar x-ray emission is potentially significant enough to increase the ionization fraction in 
nearby dense gas, potentially affecting the local gas chemistry. In this region, the complicated line 
structures found in both species near the protostars make accurate column density measurements 
difficult, however, with resulting large uncertainties. 



In their models of collapsing prestellar cores, lAikawa et al.l (|2005l ) found that NH3 can be 
enhanced relative to N2H + in core centers at high densities (n = 3 x 10 5 — 3 x 10 6 cm -3 ) due to 
dissociative recombination reactions of N2H + and e~ to form NH and N. NH then reacts with 
and H2 to form NH3. The N2H + recombination reaction is dominant only where CO is depleted 
(if CO is abundant, N2H + is destroyed mainly through proton transfer to CO), and occurs when 
the abundance ratio of CO to electrons, n( CO) / n P < 10 3 . When a collapsing core reaches higher 



central densities (n ~ 10 7 cm 3 ), however, lAikawa et al.l predict that the A^(NH3) will begin to 



decrease towards the core centre due to depletion, while iV(N2H + ) continues to increase, leading to 
a higher fractional N2H + abundance relative to NH3 at later times. This prediction is in agreement 
with X(NH3)/X(N2H + ) results since B2 is more evolved than Bl, having already formed protostars, 
but seems inconsistent with the common slope we find for the decrease in X(NH3) and X(N2H + ) 
with A r (H2) in both Cores. Detailed modelling of the physical and chemical structure of the (non- 
spherical and clumpy) Cores is beyond the scope of this study, but would help to constrain the 
central d ensity and abu ndance structure as a function of Core radius needed for a direct comparison 
with the Aikawa et al. models. 



5.4. Are N 2 H+ 1-0 and NH 3 (1,1) tracing the Oph B Core interior? 

The hyperfine structure of the NH3 (1,1) inversion transition allows the calculation of volume 
density n ex as in Equation [IJ with the opacity of a typical hyperfine component r hf = 0.233r. The 
resulting mean density n eX)NH3 = 7.8 x 10 3 cm and Ti eX} NtL 3 — 1-7 x 10 cm for Oph Bl and 
B2, respectively (note that these values are slightly less than reported in Paper I due to using the 
larger 18" beam FWHM). These volume densities are a factor > 20 less than those calculated in 
^4.1.21 from N2H + 1-0 emission (~ 26 in Bl and ~ 18 in B2). T here is some quest ion whether 



the NH3 volume densities are accurate, as recent studies (Paper I; iFoster et al.ll2009l ) have found 
NH3-derived n cx values which are an order of magnitude less (10 4 cm -3 versus 10 5 cm -3 and 
greater) than volume densities determined from continuum emission. One possible source of error 
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is the collisional de-excitation rate coefficient, j u i, which Foster et al.l note is reported in various 
publications with factor of ~ 10 variation. Since n cr = A u \ j j u i, the critical density is therefore 
also suspect within a factor ~ 10. 

It is clear that in Oph B NH 3 and N 2 H+ are tracing different motions in the gas. The gas 
traced by N2H + emission is significantly more quiescent, as has been observed in high density gas 
in other star forming cores. The fact that higher densities are calculated from N2H + emission 
relative to those calculated from NH3 emission is thus reasonable, and bolsters the hypothesis that 
NH3 (1,1) emission is not tracing the highest density gas. 

How well does the N 2 H+ 1-0 emission trace the highest density gas? In agreement with results 
in Paper I, we find that the clumps found in N2H + 1-0 emission also do not match well continuum 
clump locations, and significant offsets between N2H + 1-0 and continuum emission are also apparent 
in the integrated intensity maps. In particular, both the single-dish- and interferometer-integrated 
N2H + 1-0 intensity towards the continuum clump B2-MM8 peak offset to the clump, and in the 
interferometer data we see an integrated intensity minimum at the MM8 location. On larger scales, 
we find a trend of decreasing N2H + abundance with increasing H2 column density, described above, 
with a slope in Oph B2 that agrees within uncertainties with that determined for iV(NH3) versus 
iV(H 2 ) for both Oph Bl and B2 (see g5j] and £ j5"X4j) . Furthermore, the mean X(N 2 H+) for 
continuum clumps alone in Oph B is lower than the AT(N2H + ) over the whole core (see §5.2.ip . 
Together, these findings suggest that N2H + 1-0, though doing a better job than NH3 (1,1), may 
not itself trace well the coldest, densest material in Oph B. Appropriate caution must be made 
when using even N2H + 1-0 to probe dense gas in some star-forming regions. 



6. Summary 

We have presented Nobeyama Radio Telescope (18" FWHM) and Australia Telescope Compact 
Array (8" x 5" FWHM) N2H + 1-0 observations of the Ophiuchus B Core. Our main results are 
summarized below. 

1. The integrated N2H + 1-0 intensity is smooth at 18" resolution, and follows generally the 
850 /iin continuum emission in the Core as expected from studies of isolated regions, but significant 
offsets are found between continuum emission and N2H + integrated intensity peaks. N2H + clumps 
identified through the 3D clumpfind algorithm correlate well with those previously found in 
NH3 emission, but do not correlate well with continuum objects. Closer correspondence between 
continuum and N2H + clump locations is found in Oph Bl compared with B2, suggesting N2H + is 
better tracing dense gas in this sub-Core. 

2. We find small-scale structure in B2 through high resolution ATCA N2H + 1-0 data. Towards 
the continuum object B2-MM8, we find a broken ring of integrated N2H + emission suggestive of 
N2H + depletion. The line widths of the N2H + emission at high resolution are universally narrower 
by factors > 2 than seen in the single-dish data. Some small scale N2H + structures show velocity 
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gradients which, if due to rotation, are substantially larger than typically seen in star forming cores 
at lower spatial resolution, but are similar to other interferometric results. 

3. Over all H 2 column densities, line widths in Oph B2 are dominated by transonic non-thermal 
line widths, while non-thermal motions Oph Bl are subsonic. This result shows that N2H + and 
NH3 are not tracing the same material in Oph B, as NH3 line widths were shown in Paper I to 
be large and supersonic in both Cores. Attempting to determine a kinetic temperature by forcing 
equal non-thermal contributions to the NH3 and NjhU line widths produces unphysical results. 

4. We find double-peaked N2EU line profiles in Oph Bl. The blue line component is found 
to match the vlsr of the ambient cloud as traced by NH3 off-Core emission, while the vlsr of 
the red line component matches the NH3 Core gas. The v LSR of optically thin C2S emission also 
appears to match better that of the ambient cloud, suggesting that Oph Bl is accreting gas from 
its surroundings. Accretion may cause the co-located continuum clump (which appears itself to be 
gravitationally stable) to become unstable and form a low mass protostar. 

5. We find a larger NH3 fractional abundance relative to N2H + in starless Oph Bl compared 
with protostellar Oph B2. This result is in agreement with other studies. Chemical models of 
collapsing cores produce relatively high iV(NH3) at moderate densities (n ~ 10 5-6 cm" 3 ) and 
relatively high A(N2H + ) at high densities (n ~ 10 6 ~ 7 cm -3 ), in agreement with our results. 

6. We find a trend of decreasing fractional N2H + abundance with increasing H2 column 
density, as traced by 850 /im continuum emission. Together with relatively low X(N2H + ) at the 
positions of the continuum clumps from the single-dish data and the absence of N2H + emission at 
the positions of continuum clumps from the interferometer data, this behaviour calls into question 
the utility of using N2H + (and NH3) to trace the coldest and densest locations of star- forming 
regions. 
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Fi g. 1. — a) 850 //m contin uum emission in the Ophiuchus B2 core in Jybeam -1 , first mapped 
by iJohnstone et al.1 ( 2000bl ) at the JCMT. Black contours begin at O.lJybeam^ 1 and increase 
by O.lJybeam -1 . The 15" beam is shown at lower left. Grey and white circles show locations 
of ATCA observations with a primary beam diameter of 34". In both plots, stars indicate the 
positions of Class I objects ([Enoch et al.ll2008l ). b) The Ophiuchus B core in N2H + 1-0 integrated 
intensity in K kms -1 (TT) observed at Nobeyama. Oph Bl, B2 and B3 are labelled. (Oph B3 
is not prominent in integrated intensity due to the extremely narrow lines found at this location) 
Contours begin at 1.5 K kms -1 and increase by 1.5 K kms" 1 . The 18" beam is shown at lower 
left. Triangles indicate peak positions of N2H + 1-0 clumps dis cussed in A3. 1.11 and lis ted in Table 
[2] while square indicate locations of 850 ^m continuum clumps (jJorgensen et al.l l2008). 
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a (2000.0) a (2000-0) a (2000.0) 

Fig. 2. — N2H + 1-0 channel maps showing the locations and CLUMPFlND-determined FWHM in 
R.A. and Decl. of the N2H + clumps at the velocity channel closest to their fitted f LSR (where 
the velocity channel spacing is 0.05 kms -1 ). Greyscale ranges from -0.3 K to 1.9 K (Tmb- In 
locations where multiple velocity components were found, clumps are shown at the fitted v LSR most 
closely matching the CLUMPFlND-determined velocity at the peak line intensity. Grey contours 
show 850 fim continuum emission at 0.1, 0.3 and 0.5 Jybeam- 1 (15" FWHM). 
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Fig. 3. — a) Minimum angular separation in R.A. and decl. between indiv idual N?H + 1-0 clump s 
identified through clumpfind and 850 /xm continuum clumps identified by lJorgensen et al.1 ([2008). 
In both plots, the dashed circle has a radius matching the 18" FWHM of the N2H + 1-0 observations. 
On average, N2H + clumps in Oph B show a positional offset from continuum clumps similar to 
that found for clumps identified in NH3 (1,1) emission in Paper I. b) Minimum angular separation 
in R.A. and decl. between individual N2H + 1-0 clumps and NH3 (1,1) clumps identified in Paper 
I. Good correlation is found between N2H + and NH3 objects in Oph B. 
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Fig. 4. — a) N2H + 1-0 integrated intensity in K kms -1 (Tg) observed at the ATCA. White circles 
show the primary field of view (34") of the ATCA at 93 GHz for each of the three pointings included 
in this mosaic. The ~ 8" x 5" FWHM synthesized beam is shown at lower right. In all plots, black 
contours begin at 0.6 K kms -1 and increase by 0.6 K kms -1 . Light grey contours show negative 
integrated intensity at — 0.6K kms" 1 , — 1.2 K kms -1 a nd — 1.8K kms -1 . The N2H + clump B2- 
N5, the contin uum clump B2-MM8 (jMotte et al.lll998l ). and two Class I protostars Elias 32 and 
33 (jEliasj [1978J) are labelled, b) High resolution v LSn near B2-MM8. Colour scale is in kms -1 . c) 
Fitted Av near MM8. Colour scale is in kms -1 . 
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Fig. 5. — a) Single pointing N2H + 1-0 in tegrated inten s ity in the east Ophiuchus B2 core towards 
the NH 3 (1,1) object B2-A7 identified by lFriesen et al.1 (|2009l ). Colour scale in K kms -1 (T B ) ob- 
served at the ATCA. In all plots, black contours begin at 0.9 K kms^ 1 and increase by 0.9 K kms -1 . 
Light grey contours show negative integrated intensity at — 0.9K kms -1 and — 1.8K kms -1 . The 
34" diameter white circle shows the primary beam of the ATCA at 93 GHz. The ~ 8" x 5" FWHM 
synthesized beam is shown at low er right. The N2H " 1 " clump B2-N6, NH3 clump B2-A7 (Paper I) 
and continuum clump B2-MM15 (jMotte et al.lll998l ) are labelled, b) Line velocity or v LSR . Colour 
scale is in kins" 1 , c) Fitted Av at B2-N6. Colour scale is in kms -1 . 
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Fig. 6. — a) N2H + 1-0 integrated intensity in K km s -1 (Tb) observed at the ATCA in the southeast 
edge of the Ophiuchus B2 core. White circles show the primary field of view (34") of the ATCA at 
93 GHz for each of the three pointings included in this mosaic. The ~ 8" x 5" FWHM synthesized 
beam is shown at lower right. In all plots, black contours begin at 0.9 K kms -1 and increase 
by 0.3 K kms -1 . Light grey contours show negative integrated intensity at — 0.9K kms -1 and 
-1.5K kms" 1 . The N?H+ clum ps B2-N7 and B2-N8, NH 3 clump B2-A9 (Paper I) and Class 
I protostar Elias 33 (jEliasI Il978l ) are labelled, b) High resolution v LSR in south-eastern Oph B2. 
Colour scale is in kms -1 . c) Fitted Av in kms -1 . 
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Fig. 7. — N2H + 1-0 integrated intensity in K kms -1 (Tg) observed at the ATCA. Black con- 
tours begin at IK kms -1 and increase by 0.5 K kms -1 . Light grey contours show negative in- 
tegrated intensity at —IK kms -1 and — 1.5K kms -1 . White circles show the primary field of 
view (34") of the ATCA at 93 GHz for each of the three pointings included in this mosaic. Dark 
grey contours show 850 /um continuum emission, beginning at 0.3 Jy beam -1 and increasing by 
0.1 Jy beam -1 kms -1 . The N?H+ clump B2-N9 and continuum clump B2-MM16/162733-24262 
(jMotte et al.lll998l : lj0rgensen et alJbooa ) are labelled. The ~ 8" x 5" FWHM synthesized beam 



is shown at lower right. 
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Fig. 8.— Nobeyama N 2 H+ 1-0 spectra in T MB (K) towards the Bl-Nl (a), B1-N3 (b) and B1-N4 
(c) N2H + clump locations, showing results of single and double velocity component HFS fits. For 
each clump, the best single velocity component (1) and double velocity component (2) fits are shown 
below the N2H + 1-0 spectrum, as well as the residual from subtracting the double component fit 
from the spectrum. For each clump, the spectra and both fits are offset from zero for clarity. 
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Fig. 9. — a) Line velocity or f LSH in Oph B. Colour scale is in kms -1 . In all plots, contours show 
integrated N2H + 1-0 intensity, beginning at 3.0 K kms -1 and increasing by 1.5 K kms -1 . Only 
pixels where the isolated N2H + component was detected with S/N ^ 2.5 are shown. Stars indicate 
protostar positions. The 18" FWHM beam is shown at lower right, b) Fitted Av in Oph B. Colour 
scale is in kms -1 . c) N2H + column density iV(N2H + ) (cm -2 ) determined using Equation[3j Values 
shown have been divided by 10 12 . The maximum column in Oph Bl is iV(N2H + ) = 3.3 x 10 12 cm -2 , 
but the scale has been truncated to show variations in Oph B2. In parts c) and d), triangles 
indicate peak positions of N2H + clumps discussed in ^3.1.11 and listed in Table [5J d) Fractional 
abundance X(N2H + ) = iV(N2H + ) / N(H.2)- Only pixels where the 850 [im continuum flux density 
S u ^ O.lJybeam -1 are shown. The maximum fractional abundance X(N2H + ) = 4.9 x 10 -9 is 
found towards the Bl edge, but the scale has been truncated to show variations in Bl and B2. 
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a (2000.0) a (2000.0) 

Fig. 10. — a) Line velocity t> LSR (kms -1 ) across the continuum object B2-MM8. Contours begin 
at 3.8 kms -1 and increase by 0.1 kms -1 . The Class I protostar Elias 32 is shown. The diameter 
of MM8 is ~ 30", shown by the white circle, b) Line FWHM Av (kms -1 ) across MM8. Contours 
begin at 0.5 kms -1 and increase by 0.1 kms -1 . c) Line velocity t> LS R (kms -1 ) across the south- 
eastern edge of Oph B2, near the Class I protostar (shown). Contours begin at 3.5 kms -1 and 
increase by 0.2 kms -1 . d) Line FWHM At; (kms -1 ) across the south-eastern edge of B2. Contours 
begin at 0.4 kms -1 and increase by 0.2 kms -1 . 
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Fig. 11 — N 2 H+ 1-0 spectra towards B2-MM8, B2-N5, B2-N6, B2-N7 and B2-N8 observed with 
Nobeyama and the ATCA. Nobeyama spectra have been offset from by 2K, 2K, 4K, 2.5 K and 
2K for clarity. The temperature scale is in Tmb (K; Nobeyama) and Tb (K; ATCA). Spectra 
towards B2-N6, B2-N7 and B2-N8 are taken from the ATCA integrated intensity maxima nearby, 
and are offset from the clumps identified in the single-dish data by ~ 5 — 30". Note that little 
interferometer emission is observed towards the MM8 location. ATCA spectra at B2-N5, B2-N6, 
B2-N7 and B2-N8 have smaller line widths by a factor > 2 than are observed with the single-dish 
telescope. The v LSR blueshift of 0.1 kms -1 observed in the ATCA data relative to the Nobeyama 
data towards B2-N5 may be a result of negative features in the interferometer data. No offset 
in velocity was found between the ATCA and Nobeyama observations towards B2-N6, B2-N7 or 
B2-N8. 



-41 - 



a) 



B1 * 
B2 O 

N 2 H + clumps □ 
850 jj.m clumps A 



o 



o • 



o 



o 



o 



.o 



o 



o 



A o 



L 
21.8 



22.0 



22.2 22.4 
log(N(H 2 )) (cm" 



22.6 



22.8 



5.5 - 



4.5 



: b) 







o 


• 





O A 




- * 


o o 

* AA 


A 


h □ 


* * 


B1 * 

B2 ', 






N 2 H + clumps □ 
850 fj,m clumps A - 



21.8 22.0 22.2 22.4 22.6 22.. 
log(N(H 2 )) (cm" 2 ) 



13.2 



E 13.0 



12.8 



12.6 



1 2.4 I 
21 



c) 



□ o 

»««' 



B1 * 
B2 O 

N 2 H + clumps □ 
850 /Am clumps A 



22.0 



22.2 22.4 

| og( N (H 2 )) ( cm ~ 



22.6 



22.8 



-9 
-9 



-9.4 



-9.! 
-10. ( 



□ 






B1 * - 








B2 ; 


- p 


O 


N 2 H + 

* 850 yum 


clumps □ - 
clumps A I 


> 


* & XQ 

* 

o 

< 


A 

> '»« Q S> 












: d) 






A - 



21 .8 



22.0 



22.2 22.4 22.6 
log(N(H 2 )) (cm" 2 ) 



22. i 



Fig. 12. — a) Ratio of the non-thermal line width to the sound speed, cr NT / c s , versus iV(H2) 
(cm -2 , derived from 850 fim continuum data) in Oph Bl and B2. Each point represents an 18" 
pixel (matching the beam FWHM). The dashed line indicates cr NT / c s = 1. Most N 2 H+ clumps 
have subsonic non-thermal motions, as does a large fraction of the N2H + emission in both Oph Bl 
and B2. Note that the velocity resolution of the data, Av res = 0.05 kms , is significantly smaller 
than the apparent cutoff in the non-thermal to sound speed ratio at a NT / c s ~ 0.5. In all plots, also 
shown are values at the peak positions of N2 H" 1 " clumps ident i fied w ith clumpfind as described in 
§3 and 850 /um clump locations identified bv lJorgensen et al.1 (|2008l ). b) Distribution of n cx (cm -3 ) 
with N(R 2 ) in Oph Bl and B2. c) iV(N 2 H + ) (cm" 2 ) versus iV(H 2 ) in Oph Bl and B2. There is 
no clear trend in A r (N2H + ) seen in Oph Bl, but in Oph B2 iV(N 2 H + ) appears to increase with 
JV(H 2 ). d) X(N 2 H+) versus iV(H 2 ) in Oph Bl and B2. 



-42 - 




2 4 6 8 100 5 10 15 20 25 

T e< (K) 



Fig. 13. — a) Comparison of v LSR from N2H+ 1-0 (hashed) and NH3 (1,1) (non-hashed) determined 
from hyperfine structure fitting in Oph B. No significant difference is seen in the distribution of 
line-of-sight velocities between the two dense gas tracers. Only data points in Oph Bl and B2 
are included, since B3 was not detected in NH3 emission in Paper I with sufficient significance 
to perform HFS fitting, b) Comparison of cr NT / c s from N2H + 1-0 (hashed) and NH3 (1,1) (non- 
hashed) determined from hyperfine structure fitting. The thermal sound speed c s is calculated using 
Tk values determined from the ratio of NH3 (1,1) and (2,2) amplitudes in Paper I. Non-thermal 
motions in gas traced by N2H + are significantly smaller than in gas traced by NH3. c) Comparison 
of line opacities r from N2H + 1-0 (hashed) and NH3 (1,1) (non-hashed) summed over all hyperfine 
components, d) Comparison of T ex N2H + 1-0 (hashed) and NH3 (1,1) (non-hashed) determined 
from hyperfine structure fitting. The NH3 T ex values have a larger spread and higher mean value 
((T ex ) = 10 K) than N 2 H+ ((T ex ) = 7K). 
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Fig. 14.— Solid black lines show the isolated F X F -> F[F' = 01 ^ 12 component of the N 2 H+ 1-0 
emission line towards the three N2H + clumps in Oph Bl which show a double peaked line profile. 
The dashed curves show the Gaussian fit to NH3 (1,1) emission, observed wih the GBT, at the 
N2H+ clump peak location (observations described in Paper I). The dotted lines show the C2S 2-1 
emission at the N2H+ clump peak, also observed with the GBT. Note no significant C2S emission 
was seen towards Bl-Nl. The NH3 Gaussian fit peaks towards the red N2H + component, while 
the C2S emission peaks towards the blue N2H + component. Vertical lines matching each species 
spectrum show the locations of the fitted t> LSR for all species, including both velocity components 
from the two-component N2H + 1-0 fit. 
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Fig. 15.— a) iV(N 2 H + ) versus iV(NH 3 ) in Oph Bl and Oph B2. In all plots, the NH 3 data have 
been convolved to a final 18" FWHM beam, and both the NH3 and N 2 H + data were regridded to 
18" pixels. Only pixels where both NH3 and N 2 H + column densities were calculated are plotted. 
While the mean iV(NH3) peak at similar values for both Bl and B2, a distribution to higher column 
values is seen in Bl. On average, larger iV(N 2 H + ) values are found in B2. The combination of the 
spread to higher iV(NH 3 ) and the lower iV(N 2 H+) in lead to a larger 7V(NH 3 ) /iV(N 2 H + ) ratio in 
Bl than in B2. b) The ratio of NH 3 and N 2 H+ column densities, iV(NH 3 ) /iV(N 2 H + ) plotted as 
a function of iV(H 2 ) in Oph Bl and Oph B2. Note that significantly fewer pixels in Bl also have 
continuum emission above our threshold flux value and are thus not plotted. 
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Table 1. ATCA Targets in Oph B2 



R.A. 


decl. 


rms 


Mosaic 




J2000 


J2000 


Jy beam -1 






16:27:24.6 


-24:26:55 


0.13 


N 


B2-A5 


16:27:29.3 


-24:27:25 


0.04 


Y 


B2-A6 


16:27:31.3 


-24:27:46 


0.05 


Y 


B2-A9 


16:27:32.6 


-24:26:56 


0.12 


N 


B2-A7 


16:27:34.7 


-24:26:13 


0.10 


Y 


B2-A8 



Associated Objects 



B2-MM16 b , 162733-24262 c 



Friesen et all (120091 ) 



Wotte et all (jl998l ) 



Jorgensen et all ([2003) 



Wad (11973); IVrba et al.1 $1973 ) 
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Table 2. N 2 H+ 1-0 clumpfind peaks and parameters in Oph B 



ID 


RA 


decl. 


FWHM (x x y) 


^A,iso 


A, main 


S/N iso 




J2000 


J2000 


(AU) 


(K) 


(K) 




Bl-Nl 


16 27 13.1 


-24 29 40.3 


4500 


X 


4300 


1.45 


2.58 


7.3 


B1-N2 


16 27 14.2 


-24 29 15.3 


5600 x 4400 


1.53 


2.72 


9.4 


B1-N3 


16 27 15.3 


-24 30 20.5 


4400 


X 


4800 


2.28 


3.01 


10.9 


B1-N4 


16 27 16.4 


-24 30 30.5 


5900 


X 


6600 


1.53 


2.63 


7.2 


B12-N1 


16 27 13.5 


-24 28 10.1 


6400 


X 


6000 


1.17 


2.00 


7.3 


B12-N2 


16 27 15.3 


-24 27 50.1 


6000 


X 


6000 


1.36 


2.23 


9.4 


B12-N3 


16 27 19.0 


-24 27 20.0 


5600 


X 


4900 


1.05 


1.85 


7.0 


B2-N1 


16 27 21.6 


-24 27 45.1 


4800 


X 


5300 


0.99 


1.86 


7.3 


B2-N2 


16 27 24.9 


-24 27 15.0 


4700 


X 


6000 


1.75 


2.61 


12.6 


B2-N3 


16 27 28.2 


-24 27 40.1 


5900 


X 


3800 


1.33 


2.60 


11.6 


B2-N4 


16 27 28.6 


-24 26 40.0 


6500 


X 


6200 


0.95 


3.07 


7.1 


B2-N5 


16 27 29.7 


-24 27 05.0 


4800 


X 


6000 


1.84 


3.46 


14.5 


B2-N6 


16 27 32.3 


-24 27 00.0 


4500 


X 


6900 


2.42 


3.43 


18.5 


B2-N7 


16 27 33.0 


-24 27 40.1 


4600 


X 


3500 


1.22 


2.45 


7.6 


B2-N8 


16 27 33.4 


-24 27 45.1 


4800 


X 


3500 


1.53 


2.36 


9.6 


B2-N9 


16 27 34.5 


-24 26 04.9 


7400 


X 


5100 


1.11 


2.48 


7.6 


B3-N1 


16 27 22.0 


-24 24 54.7 


10000 


X 


7900 


1.26 


2.27 


8.8 
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Table 3. N 2 H+ 1-0 Line Characteristics in Oph Bl, B2 and B3 



Core 


Value 


Mean 


rms 


Min 


Max 


Oph Bl 


■u LSR (kms -1 ) 


3.96 


0.15 


3.43 


4.13 




At; (kms -1 ) 


0.42 


0.08 


0.21 


0.66 




r 


3.1 


2.6 


0.2 


13 




T ex (K) 


6.2 


2.0 


2.1 


12.1 


Oph B2 


v hSR (kms -1 ) 


4.05 


0.20 


3.60 


4.44 




At; (kms -1 ) 


0.68 


0.21 


0.30 


1.30 




r 


2.0 


1.1 


0.2 


5.3 




T ex (K) 


7.2 


2.8 


1.9 


17.8 


Oph B3 


■u LSR (kms -1 ) 


3.06 


0.14 


2.93 


3.44 




Av (kms -1 ) 


0.27 


0.04 


0.20 


0.35 




r 


2.3 


2.2 


0.2 


7.6 




T ex (K) 


4.5 


2.8 


1.2 


11.7 
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Table 4. Derived Parameters at N2H + 1-0 Peak Locations 



ID 


V LSR 


Av 


r 


T 

ex 


c NT / c s 


iV(N 2 H+) 


iV(H 2 ) 


X(N 2 H+) 






kms -1 


kms" 1 




K 




10 12 cm" 2 


10 22 cm" 2 


10 -io 


10 5 cm~ 3 


Bl-Nl 


3.95 (1) 


0.42(2) 


6(2) 


5.1(1) 


0.71(4) 


8(2) 


1.6(3) 


5(1) 


0.8(1) 


(1) 


3.96(1) 


0.35(2) 


0.1(1) 


6.8(1) 


0.60 










(2) 


3.29(5) 


0.4(1) 


1(1) 


3.9(1) 


0.65 










B1-N2 


3.99(1) 


0.35(1) 


3(1) 


6.4(1) 


0.58(2) 


5(2) 






1.7(2) 


B1-N3 


3.97(1) 


0.50(2) 


5(1) 


5.7(1) 


0.90(3) 


10(2) 


1.6(3) 


6(1) 


1.2(1) 


(1) 


4.01(1) 


0.35(1) 


5.4(3) 


6.9(1) 


0.58(1) 


8.4(5) 


1.6(3) 


5(1) 


2.1(1) 


(2) 


3.58(1) 


0.2(2) 
















B1-N4 


3.88(1) 


0.53(3) 


9(2) 


4.5(1) 


1.00(5) 


12(3) 


0.8(2) 


17(3) 


0.6(1) 


(1) 


4.01(1) 


0.32(1) 


3.5(2) 


6.9(1) 


0.53(1) 


5.0(3) 


0.8(2) 


7(1) 


2.1(1) 


(2) 


3.62(1) 


0.19(8) 
















B12-N1 


3.92(1) 


0.40(2) 


2(1) 


7.1(1) 


0.67(4) 


3(3) 






2.4(2) 


B12-N2 


4.10(1) 


0.40(2) 


2(1) 


7.0(1) 


0.67(3) 


4(3) 






2.2(2) 


B12-N3 


4.03(1) 


0.55(2) 


3(1) 


5.5(1) 


0.96(4) 


5(2) 


0.6(1) 


9(2) 


1.2(1) 


B2-N1 


3.83(1) 


0.55(3) 


0.4(6) 


14.1(1) 


1.02(5) 


4(11) 








B2-N2 


3.73(1) 


0.54(1) 


4.8(7) 


5.6(1) 


1.00(3) 


10(2) 


4.0(8) 


2.5(5) 


1.2(1) 


B2-N3 


4.32(1) 


0.47(1) 


1.4(8) 


9.0(1) 


0.74(*) 


5(3) 


2.1(4) 


2.6(5) 


3.7(*) 


B2-N4 


4.11(1) 


0.58(2) 


2.7(6) 


7.0(1) 


1.05(3) 


8(2) 


3.1(6) 


2.6(5) 


2.2(2) 


B2-N5 


4.25(1) 


0.50(1) 


3.1(6) 


7.5(1) 


0.89(2) 


9(2) 


3.7(7) 


2.5(5) 


2.7(2) 


B2-N6 


4.44(1) 


0.39(1) 


5.2(6) 


6.4(1) 


0.68(1) 


9(1) 


2.7(5) 


3.4(7) 


1.6(1) 


B2-N7 


4.10(1) 


0.70(2) 


1.6(7) 


7.7(1) 


1.34(5) 


7(3) 


1.3(3) 


6(1) 


3.4(4) 


B2-N8 


4.09(1) 


0.65(2) 


1.6(7) 


7.7(1) 


1.16(4) 


6(3) 


0.8(2) 


9(2) 


2.8(3) 


B2-N9 


4.03(1) 


0.61(2) 


2.9(8) 


6.3(1) 


1.13(4) 


8(2) 


2.9(6) 


2.9(6) 


1.7(1) 


B3-N1 


2.99(1) 


0.21(1) 


3(1) 


6.3(1) 


0.20(2) 


2(1) 






1(2) 



Note. - Uncertainties are given in brackets beside values and show the uncertainty in the last digit 
of the value. The iV(H 2 ) and X(N 2 H + ) uncertainties include the ~ 20% uncertainty in the submillimeter 
continuum flux values only; uncertainties in Td and k v are not taken into account. 
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Table 5. Derived Column Densities, Fractional Abundances and Non-thermal Line Widths in 

Oph Bl, B2 and B3 



Core 


Value 




Mean 


rms 


Min 


Max 


Oph Bl 


iV(N 2 H+) (xl0 12 cnr 




4.5 


2.7 


2.1 


15 




X(N 2 H+) (xl(T 10 ) 




15 


14 


3.7 


54 




n ex (xlO 5 cm -3 ) 




2.0 


2.0 


0.4 


8.9 




o"nt / c s ( kms -1 ) 




0.71 


0.16 


0.38 


1.16 


Oph B2 


iV(N 2 H + ) (xl0 12 cm" 


- 2 ) 


5.9 


1.9 


1.0 


9.8 




X(N 2 H+) (xl0~ 10 ) 




5.6 


5.4 


1.6 


29 




n ex (xlO 5 cm -3 ) 




3.1 


2.9 


1.0 


14 




a NT / c s (kms -1 ) 




1.26 


0.40 


0.64 


2.3 


Oph B3 a 


iV(N 2 H+) (x 10 12 cm" 




1.9 


0.4 


1.2 


2.9 




n ex (xlO 5 cm -3 ) 




2.0 


1.5 


0.2 


5.1 




a NT / c s (kms -1 ) 




0.39 


0.13 


0.14 


0.65 



a The 850 fim flux at Oph B3 is less than the required threshold 
of S u > 0.1 Jy beam" 1 (§4.1.2); therefore we do not comment on 
the fractional N 2 H + abundance. 
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Table 6. Derived Parameters for 850 /im clumps and Class I protostars 



ID 




Av 


T 


T 

-*■ ex 


ff NT / C s 


iV(N 2 H+) 


N(R 2 ) 


X(N 2 H+) 


Hex 




kms -1 


kms -1 




K 




10 12 cm" 2 


10 22 cm" 2 


10 -io 


10 5 cm -3 


850 /im 


162712-24290 


3.66(1) 


1.36(1) 






2.52(1) 










162713-24295 


3.86(1) 


0.52(2) 


4(1) 


5.6(1) 


0.90(5) 


8(3) 


2.1(4) 


3.7(7) 


1.1(1) 


162715-24303 


3.93(1) 


0.50(2) 


6(1) 


5.5(1) 


0.92(4) 


10(2) 


2.0(4) 


5(1) 


1.1(1) 


162725-24273 


3.72(1) 


0.61(2) 


4.7(8) 


5.3(1) 


1.16(4) 


10(2) 


4.1(8) 


2.4(5) 


1.0(1) 


162728-24271 


4.18(1) 


0.61(2) 


2.4(6) 


7.9(1) 


1.15(3) 


9(2) 


6(1) 


1.6(3) 


3.4(4) 


162729-24274 


4.36(1) 


0.55(1) 






0.94(*) 










162730-24264 


4.16(1) 


0.57(2) 


1.6(7) 


9.6(1) 


0.94(*) 


8(3) 


2.6(5) 


3.2(6) 




162733-24262 


4.23(1) 


0.73(2) 


2.8(7) 


6.3(1) 


1.35(4) 


10(2) 


4.2(8) 


2.3(5) 


1.7(1) 


Class I 


Oph-emb 5 


3.82(1) 


0.58(3) 


1(1) 


8.3(1) 


1.07(5) 


4(5) 






3.9(5) 


Oph-emb 11 




















Oph-emb 19 


4.31(1) 


0.52(1) 


1.9(7) 


8.4(1) 


0.87(*) 


7(3) 


2.3(5) 


3.1(6) 


3.3(*) 


Oph-em 26 


3.87(1) 


1.56(1) 






2.96(1) 




2.3(5) 







Note. — 850 Atm cont i nuum clump locations taken from I Jorgensen et all (120081 ). Class I protostar locations 
taken from Enoch et all (J2009J) . Uncertainties are given in brackets beside values and show the uncertainty 
in the last digit of the value. The ./V(H 2 ) and X(N 2 H + ) uncertainties include the ~ 20% uncertainty in the 
submillimeter continuum flux values only; uncertainties in and k v are not taken into account. A (*) indicates 
that the uncertainty in the value at this position were large, even though in most cases the values themselves 
are reasonable. 
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Table 7. Mean Derived Physical Parameters for N2H + clumps, 850 fim clumps and Class I 

protostars 







Av 


T 


T 

->- ex 


°"nt / C s 


iV(N 2 H+) 


AT(H 2 ) 


X(N 2 H+) 






kms -1 


kins -1 




K 




10 12 cm" 2 


10 22 cm" 2 


10 -io 


10 5 cm" 3 


N 2 H + clumps 


3.98 


0.49 


3.2 


7.0 


0.86 


7.0 


1.9 


8.0 


1.9 


850 /im clumps 


3.99 


0.70 


3.9 


6.7 


1.27 


9.4 


3.2 


3.1 


1.7 


Class I protostars 


4.00 


0.89 


1.5* 


6.5* 


1.67 


5.8* 


2.3* 


3.1* 


3.6* 



Note. — Only three protostars (two where values labelled with *) were associated with significant N 2 H + 1-0 
emission such that physical parameters could be determined from the HFS line fitting routine. 



